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Torsionsfreie achsensymmetrische Deformation von: . 
Er _ Umdrehungskörpern und ihre Inversion. £ 
ee Von A. Timpe in Berlin-Charlottenburg. 


’ 


Der etwas verwickelte Bau der die Blastizitätslehre beherrschenden Gleichungen, bedingt durch den Tensor- 
- charakter der Grundgrößen Spannungs- und. Verzerrungszustand, läßt sich für den in der praktischen Festig- 
.  keitslehre besonders wichtigen Fall des ebenen Problems und der achsensymmetrischen Deformation. von 
u Umdrehungskörpern auf eine’erzeugende Funktion zurückführen, dieder Bipotential- 
.... gleichwng genügt. Nachdem der Fall der achsensymmetrischen Torsion in einer früheren Arbeit 
a des Verf. behandelt wurde, wird jetzt. für die torsionsfreie ach sensymmetrischeDefor- 
BE NM ation von Umdrehungskörpern der entsprechende Apparat mit der Bipotentialgleichung als Schlüssel- 
ER  gleichung gewonnen, damit zugleich die Inversionsmöglichkeit. Erkenntnis: Im Rahmen der Tensorfelder 
2. 2 spielt die Bipotentialgleichung eine ähmliche Rolle wie die Potentialgleichung im Rahmen der Vektorfelder. 


2. 2 The somewhat complicated structure. of the equation dominating the theory of elasticity, due to the tensor 


un character of the Jundamental quantities siress and strain, may be reduced by the aid of agemerating 
ee function satisfying the bipotentialeqguation, inthe practically important cases of the. two-di- 
Br er mensional problem and that of rotational symmetriy. Formerly having treated the case of torsion sym- 
ET metrical around an axis, the author now develops a corresponding amalysis for to rs ionless defor-- 


. mationofrotai ional symmetry in the case of rotational bodies, terminating in the bipotential 

-  .equation, a-companied by the possibility of inversion. Conception: To tensor fields the bipotential equation 

may have similar relations as the potential equation to vector fields. ’ 

Ce travail traite le cas tr!s important pour la pratique de la deformation de type rotationel sa 's_torsion. 

En ce cas la structure um peu compliquee des &quations de la theorie de V’dasticite, ressortant du charactere 

temsoriel des quantites fondamentales „deformation“, „tensions“, est reducible & une fonction genera- 

 trice qui satisfait VEequation bipoten tielle. Le cas de la torsion de type sımmeötrique-rotationel en. 

"=  &amt traite par V’auteur dans un travail amterieur, ich analyse correspondante est developee pour le cas de ee 
la deformation symmeötriqwe-rotatiomelle, mais sanstorsiom, et c'est aussi 
V’equation. bipotentielle qu’on trouwve pour la fonction generatrice accompagnde par la possibiliti de l’inver- 
sion. Conception: Pour les camps tensoriels l’öquation bipotentielle a le röle semblable au röle de l’&quation 


‚potentielle pour les camps vectoriels. 
0801540 c.I0:KHoe IOCTpoeHne yparHeun MaTemarnyecroä Teopnuu yupyrocTu, ÖOyCHOB- 
JeHHOe TEH30PHEIM XaparTepoM OCHOBHEIX BEJINMYUH (HanpsxeHna u nedopMmannm), MO3KeT 
ÖBITb B IIPAKTHYeCKH BAIKHBIX CIYYAAX 1N0CKOH sanauu U OCECHMMETPHYHOB MedopMaluH 2 
= mei BpameHnna cBeneHo K 00P&3 yomei dyukuun, YnoBIeTBopamımmef ÖHNOTEHIHAIBHOMY A 
 ypasuennw. BonHof us IpeicHux pa6oT aBTop HCcJTeNoBANI OCECHMMETPHYHOE Epy4eHune. 
Paccmarpusaemaü sqecp cayyali oceCHMMETP uyHoä nebopManuu Te BpameHnun 6e3 
\ - EPpy4ueHHA CBONUTCH ONATb K ÖUNOTEHINHANBHOMY YPABHeHND MIA COOTBECTBEHHOH O0PA- 
 syromei byneıumm. Bsison: ÖHNOTEHNHANBHOE YPABHEHHE HMEET IA TEHSOPHOTO TIONMM AHA- 
NOTHYHOE BHAYEHHE, KAK IIOTEHINANBHOE YPABHeHMe MIA BEKTOPHOTO HOIIM. 


2 Das Problem der achsensymmetrischen Deformation von Umdrehungskörpern zeigt, wie 
Verfasser in verschiedenen Arbeiten !) dargetan hat, bemerkenswerte Analogien zu dem des 
ebenen Spannungs- und Verzerrungszustandes 2). Dabei wurde von uns die Methode der 
erzeugenden Funktionen als besonders zweckmäßig und fruchtbar weiter ausge- 

"staltet, unter Betonung der Wahl von geeigneten, der Form der zu untersuchenden Körper an- 
gepaßten Koordinaten (Zylinderkoordinaten im Falle von Zylinderproblemen usw... .). Beim 
ebenen Problem führt, wie zuerst J. H.Miche 113) gezeigt hat, das Prinz ip-der In, 
version, gestützt auf ebene Polarkoordinaten, auf einfachstem Wege zu neuen Lösungen 
von Spannungsaufgaben. Ein analoges Verfahren wurde 1931 vom Verfasser auf der Grund- 
lage räumlicher Invers ion (mit räumlichen Polarkoordinaten) für achsensymme- 


trische Torsionsprobleme entwickelt‘). 


In vorliegendem Beitrag soll die analoge Aufgabe für die to rsionsfreie achsen- 
“symmetrische Deformation von Umdrehungskörpern behandelt werden. Als Werkzeug dazu 


ı) Math. Ann. 71 (1911), $. 480-509. — Z. angew. Math. Mech. 4 (1924), S. 361— 376. — %. angew. Math. Mech. 11 


(1931),. 8. 8—15. : 
2) vgl. A. Timpe, 2. Math. Phys. 52 (1905). 


®) London Math. Soc. Proc. 34 (1902), 8.134. 
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162° Wimpe, Torionufrio achsensymmetlsche Doformation von Umdrehungskörpern a, tdi 


dient eine erzeugende Funktion F, die der Air yschen Spannungsfunktion im Falle ds 5 
ebenen Problems analog ist; letztere genügt bekanntlich der ebenen, das hier verwendete RN 
der räumlichen Bipotentialgleichung V?V/ ’F = 0. : r Ba 


 Grundgleichungen. 
Ein achsensymmetrisch geformter Körper, der kei- 
nen Massenkräften unterliegt, erfahre eine torsionsfreie _ 
achsensymmetrische Beanspruchung. Wir führen in bezug 
auf seine Achse und einen beliebigen ihrer Punkte als 
O-Punkt für den Körperpunkt P räumliche Polarkoordi- r 
naten 7,x,d ein: 
r = OP 5 & a 


y 


72 = Winkel des Strahls OP, gemessen von 
der positiven Richtung der Körperachse 
aus; re 

© = Winkel zwischen der durch P gehenden 


und einer festen Achsenebene. 


Die zugehörigen der Deformationentsprechenden Ver- 
schiebungskomponenten seien %,,%,,%g. Dadie Deformation 
torsionsfreiachsensymmetrisch erfolgt, ist allgemein us—=0, 


— = 0, und die Formänderung wird daher) durch folgende Parameter gekennzeichnet: 
Ba Er Ne Us % F 
Dehnungen: TER IT WIR az eotgy +; 
2° (u 1 du er 

Winkelverzerrungen: RER m=tz ( — E - $ | 
ehe Dehnung: ee | 3 

ubische Dehnung: Sr mRn: (r? u,sin x +7, rmein2) EEE (2) 

RS EN } RZ du 
Drehungskomponenten: 26, = 26, =0; 2 6e = = E (ra,) — 32 a (3). 


Die Komponenten des Spannungszustandes o,, 0,, 0g; Tygs Tgr, Tr, Sind mit den Formände- 
rungskomponenten durch die Beziehungen 


,—=AA+2UE; Te UYye; 
,=A4-+2us; RT ENST R hc (4) 
o=AA+2UE; mi) 
(m +1) (m —2) ’ 2(m +1) ’ 
wo E = Elastizitätsmodul, m = Querverkürzungsziffer, 


verknüpft, so daß hier 7,g = Tgr =. 


. , VondendreiGleichgewichtsbedingungen für das Volumelement 
ist eine wegen der Achsensymmetrie der Deformation identisch erfüllt. Die verbleibenden lauten: 


00, 1,97 1 
=; Ar rs + 6%, (2 2 Al 0, “en Og + 7,,cotg x) = 0; | 

RE U: 5). 
oT Tao. 1 
Sr uw: 2% A (0, — og) cotgXx +37, =0 | 


Setzt man in (5) die Beziehungen (4) und (1) ein, so erhält man zwei Gleichungen, die man 
mit ö, = ö in folgender Form schreiben kann: 


a 1 d BN 2 
or A +2u r?sin y °y aw*rsin 2), | 
Eis: R o ae (6) 


0% A+2u sing dr rn j 


*) Vgl. z.B. A. H. Love, Lehrbuch der Rlastizität, deutsche Ausgabe von A. Timpe, Leipzig 1907. 


_ Man kann diese Gleichungen auch u unter Beachtung. des. eitonietlen Charakters BP Rune 
fe nden ‚Größen - ‚aus den für rechtwinklige Koordinaten ee ade a der 
. Elasti ent wel. Love, Lehrb. d. Elastizität S. 158 BI 


ara a 


=’ [77 4 au url Ge Gy a = : 2 


m Be örsin y= w so erscheint die Ver ? 
. Knüpfung d der Funktionen p und y in der einfachen Form ER EN us < > Ba 
ot ae, ee Ka 


In BEN es One me ON 0 "siny ar 3 NE 2 


“‘ 


en gewinnen. ‚Schreibt man 2 = 9, 


Ra 


IE Re Dutch: Gielchsetzen der beiden hieraus zu ‚gewinnenden Nerte von a 


2y { ee 
ergibt sich die 
. Beziehung: RN 0% Eee. 


Po ER; Be a 5 ER = RR "= & (s = > 38 = 5 x Bu : = iR 
ae er Bear (r Pe Ft 5% sin x 3) =0 ee a 


"entsprechend der Bekannich Tatsache, daß die. kubische Dekntng A der Potentialgleichung 


en Andererseits erhält man durch Gleichsetzen der beiden aus. (6) zu REN Werte 
SAP > Rz j z r . \ 5 F 3 . 
son ar y die Beziehung: Er ee 5 | ER 


0) 1 20% ee, Re 
en Ka 2 a Se re Be 


Le a h. die für eine 1e achsensymmetrische Strömung charakteristische Gleichung. Die Verknüpfungs- 2 


23 


bedingungen (6) laufen also darauf hinaus, daß P=Aundy= nn "or sin x als Po-. 


tentialfunktion und Strömungsfunktion der achsensymmetrischen SrönunE einer idealen F lüs- 
sigkeit gedeutet werden können). . 
«  Zweiassoziierte Funktionen 9, ykann man aus einer einzigen erzeugenden Funk- 


hr 


TEE FR EN 
Ense En EAER 


E- SEEER .n 99, einer Potentialfunktion, herleiten durch den Ansatz: 

2 FR a Nee 99 | | 
f Pa per FASINS ee ee | 
= ; N 5 N I 
=» In der Tat kann man aus dem Umstand, daß @, der Potentialgleichung genügt, folgern, x 
n daß die erste der Assoziationsbedingungen (6’’) befriedigt wird. Die zweite dieser Bedingungen. ne 


’ wird aus dem Ansatz (9) als Identität gefolgert. Aus dem Bestehen der Assoziationsbedingungen 
folgt natürlich wieder ohne weiteres, daß o und y den Differentialgleichungen (7),.(8) genügen. 


\ 


Sonderlösung. 
Die Verschiebung werde als Gradient einer Funktion f(r, y) angenommen: ER 
| af 1 of w ! 
e Te 


- EB diesem. Fall ergibt. sich gemäß (a) 


02 f 2a | 
5 ee 3 


DD 

Qu 

I 
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für die assoziierte Funktion 


ee at m ,)|= Pr’ 10 

== = ee LO 
AR r?sin X AyzAG +2; er ver En 
Wegen y= 0 wird aber nach den Assoziationsbedingungen 


(6°) E-n- 0; lo o=A=V?f= const. 


5) Ein analoger Tatbestand ergab sich bei der in Zylinderkoordinaten durchgeführten Untersuchung in meiner S.1. 


angeführten Arbeit vom Jahre 1924. Auf die mathematische Rolle der den Cauchy-Riemannschen Formeln im achsen- 
symmetrischen zweidimensionalen Problem entsprechenden Formeln (6) wird noch ZUBROEZUEGENDEBUIE sein, 


11* 


TERTART EC)” 


TE ei 


_ Anteilen, einem dilatationsfreien und einem allge 


Assoziationsbedingungen (6) genügt wird, Be BE ; er 2 
teil: / ’ | £ im; , + 7 
SETBEBE Ante I Fa 2 ar Er 
waren y=— Er ER ee 

. Die zugehörigen Wertbektandteile von 5 und 26 sind gemäß (2) und 6) = re % nn 


. | Ze + n i 
Lo 
et % + 
WERE LEN Rn Van FR A 
a 2 ) D 72.7, ® . 
A * 
j ( we P iR‘ er er 
- Ri , ER 


andererseits 20, rsiny=—car'sin y' =} eine Strömungsfunktion. 


= FI ken an ung ER 
een Lösı mittels. 


Multiplikatoren dieser Anteile sind so zu wählen, daß an die Größen 4 une 25 


| = lerr)- are Kerken a) 


he er = Ir (r + taz "37m % 7% En BR nr. 


o 
FAzEM V?F); 3 : 
a Em! k ß) f - 2F — er 


a z 
Wenn P?F eine Potentialfunktion, so ist, wie Gl. (9) zeigt, e= (r 7?F) eine Potentialfunktion, 
| 2 . 


Zweiter'Anteil: 


RE Bar ir a 
% U, SINX = ——-6, EP Sr Pr 0, dr 4) 


Die zugehörigen Wertbestandteile von A und 2 & sind gemäß (2) und (3) gegeben durch: 
AR 


MECTLIEPERAERRTE nj... me 
ie \ar\ 0x 7% öylsing 9x 32 = 2 


[ * . . ” “ . = 
Wenn V?F eine Potentialfunktion, so ist 2&, rsny=@'rsiny- 


” . = 


eine Strö- 


°% 
mungsfunktion, während A, = 0 natürlich als Potentialfunktion gelten kann. 
In den durch Kopplung der beiden Anteile entstehenden Ausdrücken 


een | 2 
AA „een, ; 07 
oV?F 

9% 
sind die Konstanten c,, c, so zu verknüpfen, daß den Grundgleichungen (6) genügt wird. Das 

m 


liefert die Bedingung: 2= BT 
i Br 


20 rsiny=(—G4+6)'rsiny: 


_ 


Das Ergebnis läßt sich also dahin zusammenfassen: 


Im Falle, der torsionsfreien achsensymmetrischen 
Deformation eines Umdrehungskörpers stelle man an 
die Spitze eine erzeugende FunktionF, die der Bipoten- 
tialgleichung,. dh, der Differentialgleiehung 


x Par EN N or ee EN. (12) 
genügt. 


ee 2 ER ET N RE u 2: = | 
{ { er : 2 ee SÜD 3 SUR 5 
er E 2=v TOT m—2r dr Ep ae = 


> EN £ er. >> TR ’ T = A Mg a 
ige Verse iebungskomponenten sind dann: 
€ Hu r - i * _oF je "m e 1 bi d fr “ 4 27 i u 
u=rV? F— 2 — 4 —— et we e% 
A SEHEN 3 rsinx 0% GR ol’ 


ee Randwertproblem. | 
"Während die von A. E. H. Love °) angegebene Behandlung der achsensymmetrischen 


Verzerrung in einem Umdrehungskörper, in Zylinderkoordinaten formuliert und gleichfalls auf Be. > 


die Schlüsselgleichung V?VY?F = 0 führend, und eb iei ül lur 
n DEI NAne hung A » enso die in unserm früheren Aufsatz durch- 
ee Be erster Linie auf Körper abgestellt ist, die von Kreiszylinderflächen, den zu- 
gehörigen Meridianebenen und den zur Achse senkrechten Ebenen begrenzt sind (vgl. auch 
ra 


Föppl, Drang und Zwang II, S. 207, 208), haben wir mit (12), (12°), wie wir im letzten Ab- 


Be en werden, einen Lösungsschlüssel von allgemeinerer Bedeutung gewonnen. 


In besonderem Sinne zunächst eignet er sich für R j i 
: € andwertprobleme, bei denen die Ver- 
schiebungskomponenten %,, w, oder die gemäß (4) durch sie ausgedrückten Spannungskompo- 


_ nenten auf. zwei Kugelflächen r=r,und r=r, und i ä i 
j j = zwei Kegelflächen x = yı und X = x. VoT- 
_ geschrieben sind. Man setzt dann eine Lösung der Differentialgleichung a), eo bias? 


monische Funktion, an in der Form ?) Er 

: SB mot bg u ae (13), 
wo 9, 9* zwei harmonische Funktionen bedeuten. Für * hat hi Tenete Bier 
mentar-Lösungstypen zur Verfügung: ER ag Fe a hier geeignete Ele- 


' DE 1 
Py,=r":P„(cosx); Y, = ri Pn (cosy); 


1 
V,;,= r”  Qm (COS AB V,= mr Om (05%); 


ER: ‘ se 
hierin bedeuten die ‚,P} (cos X)“ Kugelflächenfunktionen 1. Art, die für 0 <y =” durchweg 
regulär, die „‚Qm (cos x)“ Kugelflächenfunktionen 2. Art, de fry=0 undy=n divergieren, 
im übigen regulär sind. 


Die beiden Lösungen mit dem Faktor divergieren für r—= 0, können also nur 


ym+1 


dann herangezogen werden, wenn es sich um Hohlkugelprobleme handelt, wo r — 0 nicht dem 
der Formänderung unterworfenen Körper angehört. 

Entsprechend können die beiden Lösungen mit dem Faktor Q, (cos y) nur dann eingesetzt 
werden, wenn der Körper durch Kegelflächen begrenzt ist, die die Pole aussondern. Die Heran- 
ziehung aller vier Lösungen ist also bei einem Drehkörper möglich, der von zwei 
Kugelflächen und zwei Kegelflächen begrenzt ist, 

Setzt man nun in (13) für 9, @* Reihenentwicklungen mit den gemäß (14) gebildeten, mit 
unbestimmten Koeffizienten multiplizierten Gliedern ein, so entsteht eine erzeugende Funktion F 
mit — allgemein zu reden — acht ‚‚Freiheitsgraden“. Sind bei dem oben genannten, von 
zwei Kugelflächen und zwei Kegelflächen begrenzten Körper beispielsweise die beiden Ver- 
schiebungskomponenten vorgeschrieben, so kann man die sie ausdrückenden Funk- 
tionen in Reihen nach Potenzen von ? bzw. nach Kugelflächenfunktionen von cos % entwickeln 
und durch Gegenüberstellung der aus F gemäß (12) hervorgehenden Werte von 4, 4, durch 
„Koeffizientenvergleich“ die allgemeinen Ausdrücke für %,, u, ermitteln. Die Glieder zweiten 
und niederen Ranges erfordern dabei noch eine ergänzende Untersuchung wie im Falle der 

Air y schen Funktion. Ähnlich gestaltet sich die Lösungsexekutive, wenndie Spannungss- 


_ komponenten, die sich ja durch die Ableitungen der Verschiebungen ausdrücken, auf 


den vier Teilen der Begrenzung des Körpers vorgeschrieben sind. Der Durchführung an einem 
Beispiel von besonderer Bedeutung widmet sich einer meiner Mitarbeiter. Es ergeben sich hier 
auch — wie im Falle des vom Verf. behandelten?) Kreisringproblems — mehrere Typen von 
Selbstspannungszustän den, die mit den Mehrwertigkeiten in den Verschiebungs- 
komponenten bei den niederen @liedern zusammenhängen. Künstlich könnte man sie in dem 
Körper der oben angegebenen Begrenzung erzeugen durch Eintreiben eines von zwei Ebenen 
O-const. gebildeten Kugelkeils bzw. durch Eintreiben von einem oberen und einem unteren 


Polkegel. 


ee Se 
) a. a. O., 8.317. 
°). Die Allgemeingültigkeit dieses Ansatzes wurde in der italienischen Literatur um die J ahrhundertwende dargetan. 
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Inversion torsionsfreier achsensymmetrischer Verzerrungszustände. BZ 
Beim ebenen Spannungsproblem hat die wohl zuerst von F. Levi Civita®) fest- 


gestellte Tatsache, daß die Differentialgleichung 7?7?F=0 gegenüber der Substitution 


ee y=- EEE = (mit 0®= x? + 2) invariant ist, mancherlei interessante Feststel- 
7 2? 2’ 


üglich der Inversion von Spannungsverteilungen ermöglicht, die sich ja aus einer 
ER ne genügenden Air yschen ‚Funktion herleiten lassen °); Br 
die Feststellung, daß bei der Inversion Spannungstrajektorien wieder in ee torien 
übergehen, ein gleichförmig gespannter, z. B. spannungsfreier Rand wieder in einen glei ae 
gespannten Rand. Im Falle der achsensymmetrischen Torsionszustände gibt es, wie Ve asser 
gezeigt hat!P), eine ähnliche, auf eine erzeugende Funktion gestützte Möglichkeit einer Inversion, 
so daß ein spannungsfreier Rand wieder in einen spannungsfreien Rand übergeht. i 

Für den Fall der torsionsfreien achsensymmetrischen Deformation können wir an die hier 
eingeführte erzeugende Funktion F anknüpfen, die wir in der Form 


Fer an (13) 
“ durch zwei harmonische Funktionen 91, 9, ausdrücken. Bekanntlich geht aus einer harmonischen 


i 1 ’ . . . 
Funktion @ (r, %,@) durch die Substitution r’ = Pe a 0 = © wiederum eine harmonische 


1 1 > f 
Funktion 9° (r’, y’,@’) hervor, wenn man 9’ = FRE, = 28 e) ansetzt. Diesen Inversions- 


ansatz können wir auf jede biharmonische Funktion 1 
FERN T 0) +9%(r,% 0) ee . (13°) 


F : ne 
übertragen, wenn wir =——=r'-F ansetzen). Denn damit geht aus (13’) die 
Funktion N 

F’ ' 1 1 ! ’ & ’ ’ 
mr Mila: 2A Tales 


1 1 ‚ ’ r 1 Er ) 
ale, 0)+r (3, 2,0 


hervor. Der erste Bestandteil stellt nach der obengenannten Regel eine harmonische Funktion 
von r’, x',©’ dar; der zweite Bestandteil ist das Produkt aus r’? und der harmonischen Funktion 


9% E 4%; ), das ganze also von der Form (13°), d.h.einebiharmonische Funk- 
. 1 * ”- ” 
tionvonr’,y',0'. Im achsensymmetrischen Falle stellt alsof’=r’-F e : x) eine bezüglich 


der neuen Veränderlichen biharmonische Funktion dar. Durch die so erklärte In- 
version wird also wieder eine erzeugende Funktion für eine 
torsionsfreieachsensymmetrische Deformation gewonnen. Durch 
eine besondere Untersuchung wäre noch zu klären, wie die durch sie definierten Spannungen 
im invertierten Gebilde mit denjenigen in den entsprechenden Punkten des ursprünglichen Ge- 
bildes zusammenhängen. Von besonderem Interesse sind dabei die sog. Verzerrungskerne (Sin- 
gularitäten). 

Während die für Vektorfelder bedeutsame Potentialgleichung Invarianz zeigt in einem 
weiten Felde. von Transformationen, das im zweidimensionalen Falle durch die bekannten Be- 
ziehungen zur Funktionentheorie umrissen erscheint, ist es für die Tensorfelder der bisher be- 
handelten zweidimensionalen Probleme der Elastizitätslehre (das ebene und das achsensymme- 
trische Problem) typisch, daß sie sich auf die Bipotentialgleichung zurückführen lassen und bei 
dieser die Invarianz auf den Fall der Inversion beschränkt erscheint. Man kommt zu der Auf- 
fassung, daß im Rahmen der Tensorfelder die Bipotentialgleichung eine ähnliche Rolle spielt 
wie die Potentialgleichung im Rahmen der Vektorfelder, 


?) „Sopra una trasformazione in se stessa della equazione 4,4, =0.“ Venezia (1898) (Tip. Ferrari). 
®) 8. Fußnote 3, 


10) Z. angew. Math. Mech. 11 (1931), $S. 8—15. 
11) Michell gibt dies am Schluß seines Aufsatzes beweislos an'(vgl. Fußnote 8), 


Eingegangen: November 1947. 
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N Bemerkungen über die Methode von je Duffing. ö 
zur Integration von Differentialgleichungen. 


en Von Erwin, Pflanz in Marbach/Neckar. 


'Beider Duffing schen Löstngsmethode des Anfangswertproblems einer Diffgl. Y’(%) =f(x, y) mit | 


Ya) = Yı ergibt sich für eine Näherung Y von Y (x) (x = x, + h) eine Gleichung, in der auch Ableitunge 
von f auftreten. Diese und ähnliche Bestimmungsgleichungen Bar Näherungen Y (beliebig hoher en 
in h) von y lassen sich auch aus einer Verallgemeiner@ng der Taylorentwicklung für Funktionen einer Va« 
riablen entnehmen. Unter Voraussetzung der Monotonie einer gewissen Ableitung von f im, betreffenden 
Intervall kann man den genauen Wert y (x, + h) durch solche Näherungen Y in Schranken einschließen. 
Nach Übertragung der verallgemeinerten Taylorentwicklung auf Funktionen zweier Variabeln werden damit 


zwei Verfahren (2. und 4.Ordnung) zur Lösung partieller Differentialgleichungen 1. Ordnung hergeleitet. 


Duffings method for solving a differential equation y’ (x) = f (x, y) with the initial values y (x) = Yyı 
produces an equation for Y — an approwimation of y(a)[x = x + h]—, in which derivatives of f too, 
are appearing. This equation and other similar to it, giving approwimations Y (of an arbitrary order in h) 


‚of Y, also result from a generalization of the Taylor development for functions of one variable. Supposing _ 

the monotony of ‚% certain derivative of f in the interval concerned, the ewact value y (x, + h) can be included 
- between two limits by approximations of the kind. The generalized Taylor development is transferred to ; 

functions of two varıables, and so two methods (of the second and the fourth order) are derived for the solution 


of partial differential equations of the first order. 
La methode de Du f fing pour la solution du problöme des valeurs initiales d’une Equation differentielle 
yY’ (x) =f (8, y) avec Yy (x) = yı fournit pour une approximation Y de y (x) (© = x| + h) une öquation 


qui contient aussi des derivees de f. Cette Equation (et autres semblables) qui donne des approximations Y 


(d’un ordre arbitraire en h) est aussi accessible par une generalisation du developpement selon Taylor. 
l’on suppose la monotonie d’une ceriaine derivee de f dans l’intervalle considere, on peut enfermer en limites 
la valeur exacte y (x; + h) par de pareilles approximations Y. En adaptant le developpement selon Taylor 


ü des fonctions de deux variables, on arrive & deux methodes (du deuxieme et quatrieme ordre) pour la 
z 


‚solution des Equations differentielles du premier ordre. 


Meroa Ayhhuura pemenus nuhbepermmansHoro ypaBHeHHs Yla)=f(e, Y) c HAYANBHEM - 


ycnopnem % (x) = Yı AaerT ua mpnönmkenna Y hyHneuum %(2)(®—= x, +h) ypasHeuue, 
KOTOPoe CONep:KuT TAKIKe U NMPONSBONHBIe OT f. JITO mim IOoMNoÖHOe eMy YPaBHeHne 
(1100604 CTeNeHH OTHOCHTENEHO h) qua npnönnkenna Y hyHRUMU y Mo:Ker ÖbBITB TONyyeHo 
.TAKKE U IIPH IOMOoImMm Pasıo:keHun' PyHKIun OT OnHOoäA MepeMeHHOH B OÖ0ÖMeHHBIH pa 
Tafrnopa. MHcexona us mpennonokeHung, YTO HEKOTOpaA IPOHSBORHaAA OT f B NAHHOM 
HHTePBale HIN NOCTOAHHO BOSPacTaeT HIN HOCTOAHHO YÖBIBAeT, MO3KHO TOYHOE 3HAyeHHE 
hyarıun Y (x| + h) sakııoyuTb Meikıy AByMA npnönmerkHuuamu Y. Pacnpocrpansa pasıokeHne 
B o000meHHsä pay Taünopa Ha byHrımm OT AByX HepeMeHHEX, ABTOp ToNy4yaeT ABa 
MeTON& (BTOPOTO H YETBEPTOTO NOPAAKA) pemeHnsn nuhhepeHmnansHbIX VPABHeHHH B YACTHEIX 
HPOUSBOAHBIX HEPBOTO TOPAAKA. 


1. Einleitung. Bei seinem Verfahren der genäherten numerischen Lösung von ge- 
wöhnlichen Differentialgleichungen 1. Ordnung. 


d v 
v-,-/® y)...(la) bzw y—n=/&EYyÖ)dd.... (1b) 


t=1 


mit <= x; +h und der Anfangsbedingung y (x,) = y, ersetzt G. Duffing!) das Integral 


; h ’ 772 N . . . 
in (1b) durch den Näherungsausdruck F (4y) +hy) +2y'). Damit ergibt sich aus (1b) 
die Bestimmungsgleichung 

rent TER ZRN ENT...) 
für eine Näherung Y von y an der Stelle <=x, +h, wobei ,=Y(w)=y, ı, = 
Y/’= y/ vorausgesetzt ist. Die Auflösung von (2) nach Y ist ersichtlich praktisch einfach 


nur dann durchführbar, wenn f (x, y) linear oder quadratisch in yist. Die Entwicklungen von Y 
bzw. y an der Stelle x, nach Potenzen von h = x — x, ergeben Übereinstimmung in h bis ein- 


‘schließlich der Glieder 3. Ordnung 2), wie auch im folgenden Abschnitt 2 gezeigt wird. 


Die Genauigkeit dieses Verfahrens ist also, wenn man die Zahl der übereinstimmenden-Glieder 
bei Potenzreihenentwicklung als Maßstab für die Güte der Annäherung ansieht, gleich derje- 
nigen der Runge schen Formeln; sie ist somit größer als die des Cauchyschen Differenzen- 
verfahrens, geringer als die des Verfahrens von Kutta°). — Von P. Funk *) wurde darauf 
hingewiesen, daß „‚hier ein ähnliches Prinzip vorliegt, das auch bei einer bekannten Herleitung 
der Eulerschen Summenformel und auch beider üblichen Herleitungder Runge-Kutta- 


ı)G.Duffing, Zurnumsrischen Integration gewöhnlicher Differentialgleichungen 1..und 2. Ordnung, Forschungsarb. 
Geb. Ing.-Wesen 224 (1920), S. 29—50. 

2) Fr. A. Willers, Methoden der praktischen Analysis, S. 307. Berlin 1928. h 

s) Fr. A. Willers, Numerische Integration, S. 89—92, Sammlung Göschen, Bd. 864, Berlin 1923. E 3 

“4 P. Funk, Über Duffings Methode zur numerischen Integration von gewöhnlichen Differentialgleichungen, 


Z. angew. Math. Mech. 7 (1927), S. 410, 


iB 
. ar 


eh 


EE, 
BEN 4 

Dar 2; 
ww 


finden lassen, 
4) (& 
hen 


zweier unabhängiger Variabeln, die Duffingsche Methode zur genäherten Lösung partieller 5 
Differentialgleichungen 1. Ordnung, wobei wir hinsichtlich z Näherungsverfahren zweiter und E% 


rhalten. E EYAk Kt 

© Fear Duffingschen Methode bei gewöhnlichen Differential- 
gleichungen 1.Ordnungy=f(a,y)mit der Anfangsbedin zungy= ern: 
‚für 2=x,. Ausgangspunkt unserer Betrachtungen ist die verallgemeinerte Taylorssca 
Entwicklung ) | — a = 


a Br > A he a a 
ur r , 7 
4° F . i” Er 


Be _s0 Konya DB LA.R 3 > 

4A [® (x) — B(z,)]= & A ee? (%,) ( ) x | (2) + x a ( ) E F- 

’ „1 ul n=1 ul 2 

3 mit \ a ar 

Be & ; ». Edle) R : Fa 

Ba N a De ae = ( = PR ee — 
; und A f qyP ® ü 

NEN El N en B n(%+B+2D) = 


s 3 
Aus (3) wird, beia=0 bzw. ß=0 die gewöhnliche Taylorentwicklung, wenn man & bzw. 


[d p=1 
& wegläßt und A=1, Ar 0 setzt, 
u=1l 
e 3 RE, 
Mit B (x) = y (a), Y=f(% y), Y" = — +-—- :f=9 (%y) erhält man aus (3) für 
ee 2, Bß=1 E 0% oy 


ya yon ty ty — See HyOldar >. (da) 
und für BZ En . 
Aha) AHRgEY)=6yn +2 + EEE HyMd)aE.. (4b). 


t=1 
Vernachlässigt man in (4a), (4b) die Integralausdrücke, so ergeben sich Bestimmungsgleichungen 
für zwei Näherungen a bzw. y4® von 4 (x). Dabei sei vorausgesetzt 


ven = en — y,, yan' ht Ku er y. ER, ni y* 
Man erhält 
erm—htar] =6y+4Ahyy+hy’ ... (a) 
6 y (12) A hf(e, le +Ah2y(e, y9, = Yyı #33 y| a (5b). 


(5a) ist die Duffingsche Gl. (2). 


Darstellungen für die Abweichungen nad (2) = Y&V _ Y, zy (2) = rd y ergeben _ 


sich durch Subtraktion von (4a) bzw. (4b) von (5a) bzw. (5b) und Anwendung des Mittelwert- 
salzes: 


°) Wir verwenden hier im wesentlichen die Darstellungsweise von N. Obreschkoff, Abh. pre 
h = i 2 5 B . uß. Akad. Wiss, 
Math.-naturwiss. Kl. Nr.4, 1940 S. 1— 20, Diese Verallgemeinerung der Taylorschen Formel findet er samt a 
schon bei S. Hermite, Oeuvres III (1912), 8.438. Hermite geht bei der einen seiner beidenHerleitungen aus von der 


Henlskhne [a „(m) (a , Re BEA un a (y®—t, m) BE EC 


£: ß “ 
Mitu = BT v=(2—a) (c—b) ‚(,ß=0,1,%..,08+Pß= m) gewinnt H. hieraus wegen A) = m! eine Quadra- 


turformel für [ (2)4 8, die beid'(!) =p(L}), a = x, db = z gleichbedeutend ist mit (3). Weitere Beweise zu (3) sind 
[20 

außer beiN.Obreschkoff (8.o.) noch angegeben bei G.Kowalewski 

und L.. Tschakaloff, Ann. Univ. Sofia, Fac, Phys.-Math. 34 3 

Math. 26 [1942] 8. 11), 4 Be 


Dtsch. Math. 6 (1942), 8.349 bis 351 
8. 353—394 (letzteres Zitat nach Zentralbl. £. 
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“ 
% 


nr [2 a eo] = [en Ars" nar- 
2a 0) 
= ya +9 n) 


nannte) 4128 a He] = 


ar | 2 (4) d ht (4) (1,2) 
So Ve Fe) (& +4 h) 


== % . 
(er, 9], lea am A < 1) < 
Offensichtlich folgt aus (6a), (6b) durch Reihenentwicklung für kleine || für die Abweichungen 7: 
Reh hi 
2 = 53479 y®(z (% 1) + ...9 (7a), 27 = 29 y® (x (x 1) TE BET (7b), 

und für die Abweichung A der Y: 

- 4 

RE RE RER Li RIP: 


Aus (7a), (7b) ersieht man, daß (5a), (5b) Näherungen Y liefern, die mit y bis einschließlich der 
Glieder mit A? (bei Entwicklung nach Potenzen von h an der Ausgangsstelle x,) übereinstimmen. 

Aus (6a), (6b) folgt, daß für genügend kleine |A| n”® und n® ’) verschiedene Vor- 
zeichen besitzen, d.h. Y'*" und Y"? den genauen Wert y einschließen, wenn y»(d) in 
(2), % +h) sein Vorzeichen beibehält. Ob das gewählte |h| „genügend klein‘ ist, ist 
durch Betrachtung der Klammerausdrücke auf der linken Seite von (6a), (6b) zu entscheiden. 


Für die Abweichung Y„— y des Mittelwertes „(ren Tr y@2) der beiden 
Näherungen folgt 2 


ee erde h5 2 (8) HI Be) 
Im y-,M Re so 7 @)r5y A, ur nn) 


af (8). 
2] +... 

Y„ nähert also y etwa mit der Güte der Kuttaschen Formel an. (8) ergibt sich, wenn 
man in (7a), (7b) noch die Glieder 5. Ordnung anschreibt und dabei in den (wegen (6a), (6b) 


rechts) auftretenden Integralen im Integranden y®%Xc) an der Stelle &= x, entwickelt. 
Beispiele: a) Die Differentialgleichung laute Y"=ymit y=- y,=1für z,=0; es sei 
== 0:L, 
i 6,41 6,2 
E - (2,1) „a 
Dann folgt aus (da) Y \ 58 561 
‚Da yıd) sowie für das gewählte = 0,1 die Klammern auf der linken Seite,von (6a), (6b) 
ihre Vorzeichen durchweg beibehalten, gilt 1,105 169... < (0,1) < 1,105172... 
Weiter folgt Y„ = 1,105 170 87... 


Der genaue Wert ist y (0,1) = e”' = 1,105 170 91. 


—= 1,105172..., aus (bb) 7” = — 1,105.169. 


4 
b) Bei dem D uffingschen Beispiel ®) en Helfen erg 


ergibt sich aus (da) 7" = 0,751 676..., aus (5b) 7” = 0,751 017... 
Weiter gilt, wie man aus (6a), (6b) unschwer nachweist, 
ZUR ey iL,1)<70,781676 2% 
Mittelbildung ergibt Y,„ = 0,75134..., die Kuttasche Formel?) 0,75121...., während 


1 
der genaue Wert der Lösung y = SE Il YyeED = 0,2131... ist 


Aus (3) lassen sich [auf gleiche Weise wie (5a), (5b)] Bestimmungsgleichungen für Nähe- 
rungen Y beliebig hoher Ordnung in h herleiten; man braucht in (3) nur « + ß genügend groß 


zu nehmen. Gewisse Paare 7”, Y®* von Näherungen schließen y ein, wenn |h| genügend 


6) @&., Duffing,a.a.0.:). 
’)_ Fr. A. Willers,a.a. 0.°) 8.9. 


- Sala 
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klein ist und y*+?+D(e) im jeweiligen £- Intervall (x, x, +h) sein Vorzeichen 

wechselt. 
So ergibt sich für das eben betrachtete Beispiel a) bei «= 3, P= 2 ne 
y®® _ 1,105 170 91792... <y(0,1)= e”! < 1,105 17091823... Y®”, 


wobei y (0,1) durch Y®®, y@® in 5. Ordnung, durch den Mittelwert =, (r®” + y®®) in 


nicht _ 


sechster Ordnung angenähert wird, Dieses Verfahren der Einschließung und Mittelwertbildung 
kann gelegentlich bei Auflösung linearer Differentialgleichungen 1. Ordnung von Vorteil sein. 
Auf die Herleitung entsprechender Formeln aus (3) zur Lösung von Differentialgleichungen 
2. Ordnung und deren Vergleich mit den Ergebnissen von Duffing und Funk sei hier 
verzichtet. 

3. Übertragung des Duffingschen Verfahrens auf partielle 
02 0z i 
a =) =0. Wir geben 
zuerst die (3) entsprechende Erweiterung der gewöhnlichen Taylorschen Formel für Funk- 
tionen 2 von zwei unabhängigen Variabeln. Hierbei werde die gleiche Methode herangezogen, 
mit der sich auch die gewöhnlichen Taylorentwicklungen von Funktionen einer bzw. zweier 
Variabeln in einfachster Weise ineinander überführen lassen ®). 

Nach (3) gilt, bei anderer Bezeichnung der unabhängigen Veränderlichen, 


Differentialgleichungen erster Ordnung r(, Y, 2, 


& 1 ß 
APLS (t) — © (0)] = ZA, 2 Pe gm () +AP-R 
uh= 


&—H De u! 
mit 
BE al [ur gere+rn (7) ds 
(& +B)!o 
el a 


ET NN En ortP+D »t ; 0 1: 
pi (AH HI! SEE 
Für d()=2(, +h,t,y, +) (., HZ C— u. = Yy— y, fest) 


ergibt sich hieraus, wenn man die Ableitungen nach £ bildet und nachträglich noch ? = 1 setzt, 
die gesuchte Verallgemeinerung 


& fl 0) o\w) 
Met net= (m rn) 
EN TE era er Ri 
B u (u) vn 
BR +1) h 2 2) ß 
ZA (th) Ban HA 
® (1 al pl ( ? O\la+s+1) 
mit ROH HBFTI hs, +%, 2 2 (%, +9, Rh, Y%, +9, k,), VEIT 


Dabei bedeute in symbolischer Schreibweise 
ß) ß) (u) 7} u x 
(r +62) 2 (®, y) = h" 0 2(2, y) N a} "2 (x, y) 


0% 0 dx" aa tay 
Ua LE PA) 
a ay® öy* 


Für “ =) bzw. ß = 0 werden aus (9), unter Beachtung des im Anschluß an (3) Gesagten, die 
gewöhnlichen Taylorentwicklungen von z an den Stellen (z, y) bzw. (x, %,) 


w 


Zur numerischen Bestimmung von Punkten (z, y,2)derdurch die Anfangskurve Glz, y;2)=0, 
H (%, 9,2) = 0 gehenden Integralfläche ®) der partiellen Differentialgleichung 


Past A a WR 
05 Y, nr Fl dy En al a an (10) 
werden wir unter Heranziehung von (9) zwei Verfahren betrachten. 


°) R.Courant, Vorlesungen über Differential- und Integralrechnun 
h ” g 8, 2. Bd., 8.62. Berlin 1929 
°) Wir setzen dabei voraus, daß durch die gegebene Anfangskurve sich genau ei 5 
a vuß, u eine Integralfläche legen lasse 
die in den folgenden G@1.(12) und (16) auftretenden, zu Punkten dieser Anfangskurve gehörigen partiellen longer 


eindeutig bestimmbar sind. Weiter seien, wie auch sonst in diesem A. isti 
a em Aufsatze, alle angeschriebenen Ableitungen existierend 


Für ur re Er rs au der ar a % Y ee — I H he 2, we A Ze die wir uns gege. 
on, kann c, nach (12) als bekannt angesehen werden. Dann hat man, wenn man (1 
Anfangskurve anschreibt, zusammen mit (10) 3 Bestimmungsgleichungen 


RE 92 (w, | 
‚ herungen ( der 3 Unbekannten 2 (® y), Kae Be (Der Einfachheit halber se KEN 


NO Rd re 


2% es Näherungen keine besonderen ‚Bezeichnungen hier eingeführt), 


5 ‚2 (. ir ist zu bestimmen aus 


Man erhält Er a Be Ben 

ER ON er te, IT 
| ea, in, = 18" = WB n .. . a) 

ET & ae = Re 2% „h 


> 


Flay 92, (a 2) Er PN=0. EL a ee aa. 
£ FRE } fa hre n II. Mta=ß=2 ergibt sich aus (9), wenn man u R vernach- 


lässigt, 
netmn-olk oo +4, EIER RIEHEN 
Ba * = NO 
2 | | MY): 92m en Be 
‚+2h, PR ITE Rare =d,. 
‚mit | RE jr 2 (z,. %,) 0% (2, y,) 
RD R 6 1, 'v» % v. | | 2 v». Yo 
era 
f 0 ?2 (2%, y,) 2 02 (0, 9) I: 5 
*£ 2 h k, ey: k} TT dyRE- ... 
Schreibt man (15) für drei, auf der Anfangskurve liegende Punkte (&, %ı): ... an, so.hat 
man zusammen mit den Gleichungen 
02 02 
Plannn= 2, re 17) 
PFAFF DHF, tr HF, =0 i 
Pu +Rsg+Fy +, =0 | 
2 02 92% 
6 Gleichungen für Näherungen der : men Pe ,yen sen 
022 022 0 9y Das: 
= ad’ I an der Stelle (z, y). Be 
Nach (16) ist d, als bekannt anzusehen. 
Aus den 3 linearen Gleichungen (15) @=1,2,3) kann ma Dr 002 durch 2, 
KLEE 8 ER r BurrT 080y’ Ay? eek 


— ausdrücken und in (17) einführen. - Aus den so sich ergebenden Gleichungen ist 2, 3 


zu bestimmen. 
Bo ineneatygeklingen unter Heranziehung von (9) zeigen, daß Verfahren I bezüglich z von 


2, Zaun lc ae von 1. Ordnung, Verfahren II bezüglich z von 4. Ordnung, be- 
02 EU Bo 2 9z 
5 von rdnung, bezüglic da’ Dudy’ dgp 
a geht also bei Verfahren I die Abweichung des Bier bei sich ergebenden Nä- 
herungswertes für z vom genauen Werte mit 0°, bei Verfahren II mit 0° gegen Null bei 0 0, 
wenn man mit o etwa den Radius der kleinsten Kreisscheibe bezeichnet, welche die Stelle (x, %) 
und die Punkte (x,, Y,) in ihrem Innern oder auf ihrem Rande enthält. 


züglich 2° von 2. Ordnung ist. 


s=h, h, N h, & 


nn ” re Y.)7, ex an 2 er: ‚ = ku Fl 4, 


a Fin; Kr E 2 Zur Veranschaulichung a 


Zu bestimmen sind hierfür Punkte z des partikulären Integrals 


Für die Stelle x = 1,1, y=1,1 ergibt die genaue Lösung 
FE | Da | 92 ®2 
2 u Be: 
| 0% 


02z erg 


1,107 090 6,.. . 0,940 735 8,. . 0,221 418 15... 


3 | 32 bi 
2 | dx [027] F 
108 MIR, 0,2200... RE 
Abweichung: 
—0,9%/0 —2, 20/0 6,3900 A 
und nach Verfara U =, = 3-1, „-1,%=L1, „= 1,2) aus (15) bis (17): 
02 02 022 2% = . #z 
< °x 9y 9x2 0x 0y oy? 
BEER TER FR E SEBER AR SN ORNE GE ERHRIK EEE EEN TEEN DH u nn 
1,107 087 49...| 0,940 692... 0,231472: 7A 241I IN HR 778 0,252 699... 
Äbweichung: ı 
0,0039 0,0496 — 0,2490 — 1,6% re —2,80/, 


Ohne auf die Verhältnisse bei partiellen Differentialgleichungen 2. Ordnung näher ein- 
gehen zu wollen, sei noch bemerkt, daß z. B. zur Lösung der Differentialgleichung 


2 ?2t ) 
RA.) er Ne u . + (18) 


0% ay? > 
“bei vorgegebenen Randwerten 2 (2, Y,) sich aus (9) speziell mit «=0, ß=2 (bei Vernach- 
' lässigung des Restgliedes), also durch Anschreiben der gewöhnlichen Taylorentwicklungen 
von 2(2,9) ®=1,2...5) an der Stelle (x, y) bis einschließlich der Glieder mit den Ab- 
leitungen zweiter Ordnung, zusammen mit (18) ein System von 6 linearen Gleichungen für 


%2 (x, y) 02 (x, Y) . . 
FR aufstellen läßt. Hieraus ergibt 


Näherungen der 6 Unbekannten z (z, %), 
sich als Näherung für 2 (x, y) 


1 
gar... +2)... en RE SR Er 


In (19) bezeichnen 2, diegenauen Werte zan fünf, gleichmäßig auf der Periplierie eines Krei- 
ses mit Radius o und Mittelpunkt (x, y) angeordneten Stellen x,, %,, wobei 3 =%, „=y-+o sei. 

Berücksichtigt man in (9) (bei «= 0) auch die Glieder höherer als 2. Ordnung, so erhält 
man aus den eben erwähnten, aus (9) folgenden, 5 linearen Gleichungen die Beziehung 


Pzlmy) Men) A E 3.0? [etz (z, y) 
Im 173m, Falten ee 
+2 oz(m,y) Otlz(m,y) ur 


am | Day SR 


Mit (18) folgt aus (20) ersichtlich wieder die Näherung (19) 10), 


1°) Im speziellen Falle (& = 0) der gewöhnlichen Taylor entwieklung kommt man also hier im wesentlichen auf das 
bei Aufstellung von Differenzenausdrücken (als Annäherungen an Differentialausdrücke) übliche Verfahren. In (20) haben 
wir mit dem ersten Klammerausdruck rechts vom Gleichheitszeichen einen solchen Differenzenausdruck für 42(x, y) 
bei Heranziehung der Funktionswerte in N = 5 gleichmäßig auf dem Rande eines Kreises um (z, y) mit Radius o verteilten 
Stellen. Im Falle N = 3,4, 6 finden sich diesbezügliche Ausdrücke für Az samt Angabe von Fehlergliedern bei L. Col- 
latz, Eigenwertprobleme und ihre numerische Behandlung, 8. 286, 293, 296, Leipzig 1945. [Zusatz bei der Korrektur: 


Wie Verf. nachträglich bemerkt, sind der Gl, (20) entsprechende Beziehungen für all i 
S. A, Gerschgorin, Isw. Leningrad, Polytechn. Inst. 1927) = a 


angegeben bei Ch. { ; 
s. Nr. 6, 819841 (1934). DE ei Oh. Mikeladze, Bull, Acad. Sci, U,R,S,8, VI, 


Eingegangen am 26. Oktober 1947, 


Sag 3x = a. 
a 
Ti 1.2 me | ee 
| 2,423 366,...| 3201039... | 0,245572... 
Nach Verfahren Imty =, = 1, „=1, % = 1,2) erhält man genähert aus (13) aa): ix % 


Ye 


Runge-Kutta-Verfahren zur numerischen Integration 
von Differentialgleichungen n=ter Ordnung ne 


Baer Now Rudolf Zurmühl in Darmstadt. a 


0. Im Verallgemeinerung der meist benuizin Runge-Kutta-Formel für Differentialgleichungen. 
4. Ordnung wird ein Integrationsverfahren nach dem Runde En ita- nn a I 
..... gration von Differentialgleichungen beliebig hoher Ordnung angegeben, welches sich gegenüber dem gewöhm- 
0.00 lichen Runge-Kutta- Verfahren durch eine mit der Ördnungszahl n wachsende Genauigkeit und durch 
00. Arbeitsersparnis in Sonderfällen auszeichnet. Das .Verfahren, welches insbesondere zur Integration von 
Differentialgleichungen 4. Ordnung vorteilhaft benutzt werden kann, wird an einem Zahlenbeispiel einer 
solchen Differentialgleichung erläutert. AETSER VDE RE a na En Ne ER 
> The mostly used formula by Runge-Kutta apt for differential equations ‘of the first order is 
gemeralized ihus a method following the Runge-K BE Deine is Stated to re in one step 
. differential equations of an arbitrary order n. Computed, with the usual Runge-Kwutta method, the 
given ‘method is distinguished by an accuracy inereasing with the number n and by the possibility of saving 
computing work in certain special cases. The method especially advantageous in integrating differential 
equations of the fourth order is illustrated by a numerical example of a differential equation of that kind. 
 .. ...  Ilyrem o6oömenns hopmyssr Pyure-Kyrra, ynorpedamomeica yanıe apyrux a aude- 
.  — PeHIHANBHBIX YpaBHeHnä LePBOrO TIOpAnKa, IONyyaeTca Ha ocHoBe ıpuunmma Pynure-Kyrra 
5 . 0. MeTON HEIOCPEACTBEHHOTO HHTETPHPOBAHHA AUPePeHNHABHEIX yPaBHeHuH MO00TO TOpAanka. 
NEE: Io cpassueumo c oösunsım meronom Pyure-Kyrra oH orımyaerca 0codeunoi mpocroroh 
Re B HEKOTOPEIX YACTHBIX CHYYAAX A TAKIKe TeM, YUTO ETO TOUHOCTB BOBPACTAeT C TOPANROM 
 n AuhepeHlMasHoro YPABHeHHA. ITOT METON, OCObeHHO MHPHTONHBIÄ AI HUHTerPHPoBaHnua a 
-  AudepeHmmassHbIX YPaBHeHMÜ YETBePTOro HOPAAKA, HIMOCTPHPYETCH HA YHCHeHHOM Dpnmepe 
TaAKOTO ypaBHeHun. 


1. Einleitung. = 

| Die numerische Integration gewöhnlicher Differentialgleichungen höherer Ordnung läßt 
> - sich bekanntlich in der Weise durchführen, daß die gegebene Gleichung n-ter Ordnung als System 
' von nDifferentialgleichungen 1. Ordnung geschrieben und auf dieses System eines der Integrations- 
verfahren 1. Ordnung angewandt wird. In dieser Art ist insbesondere das-bekannte Runge- 
 Kutta- Verfahren zur Integration auch von Differentialgleichungen 2. und höherer Ordnung 
herangezogen worden. Indessen hat man namentlich für die in den Anwendungen wichtige Diffe- 

“ rentialgleichung 2. Ordnung auch nach unmittelbaren numerischen Integrationsverfahren gesucht. 

Das Runge-Kutta-Prinzip ist dabei zuerst von Nyström!) angewandt worden, welcher 
derartige Formeln der verschiedenen Genauigkeitsstufen zur unmittelbaren Integration von 
Differentialgleichungen 2. Ordnung und von Systemen solcher Gleichungen aufgestellt hat. Von. 
diesen Formeln ist besonders ein Formelsatz praktisch wichtig geworden ?), welcher der meist 
benutzten Runge-Kutta-Formel für Differentialgleichungen 1. Ordnung entspricht und 

wie diese eine Genauigkeit der Ordnung h* aufweist, wenn man mit % wie üblich den Argu- 
mentschritt bezeichnet. ‚Die Formel hat gegenüber der Runge-Kutta- Formel für Diffe- # 
rentialgleichungen 1. Ordnung den Vorzug einer etwas größeren praktischen Genauigkeit (vgl. 
z.B. Tafel 7 dieser Arbeit, Spalte 6, 7, 8) sowie einer Arbeitsersparnis in Sonderfällen. 3 

Von Runge-Kutta- Formeln zur unmittelbaren Integration von Differential- 

gleichungen höherer als der zweiten Ordnung ist bisher nur eine Formel für Gleichungen 3. Ord- 

nung angegeben ®). Für die praktisch wichtigeren Differentialgleichungen 4. Ordnung lag zunächst 

die Aufspaltung in zwei Differentialgleichungen 2. Ordnung und deren Behandlung nach den 
‘Nyströmschen Formeln nahe. Indessen kann auch für Gleichungen 4. Ordnung die Aufstel- 
lung einer unmittelbaren Integrationsformel versucht werden in der Erwartung, daß hierdurch 
eine weitere Genauigkeitssteigerung eintritt. Dies trifft nun, wie sich im folgenden zeigen wird, 
£ in der Tat zu, und zwar noch in weit höherem Maße als beim Übergang von der alten Runge- 
Kutta-Formel 1. Ordnung zur Nyströmschen Formel 2. Ordnung. Während dort nämlich 

- die Größenordnung h* der Genauigkeit erhalten bleibt, wird demgegenüber die Genauigkeits- 

stufe bi Runge-Kutta- Formeln für Differentialgleichungen höherer als der zweiten Ord- 

nung mit fortschreitender Ordnung wenigstens teilweise erhöht ?). Hierdurch aber erfahren die 
numerischen Ergebnisse u. U. eine so beträchtliche Genauigkeitssteigerung, daß nicht allein das 


ı)E. J.Nyström: Über die numerische Integration von Differentialgleichungen. Acta Soc. Sci. fennicae Bd. 50 


1925), Nr. 13. ; 
! Vgl. etwa R. Zurmühl: Zur numerischen Integration gewöhnlicher Differentialgleichungen zweiter und höherer 
Ordnung. Z. angew. Math. Mech. Bd. 20 (1940), S. 104—116, insbes. S. 109—112. Die praktische Brauchbarkeit des dort 
beschriebenen Verfahrens ist in mehrjähriger umfangreicher Anwendung auf Differentialgleichungssysteme 2. Ordnung 
am Institut für Praktische Mathematik der T.H. Darmstadt zur vollen Zufriedenheit erprobt worden, 
:) R. Zurmühl,a.a. 0.2), S. 113—115. . 
*) Dies zeigte sich bereits bei der vom Verfasser aufgestellten Formel für Differentialgleichungen 3. Ordnung als wobei 
für y eine Genauigkeit der Ordnung h° erreicht werden kann; doch wurde dieser Umstand damals vom Verfasser wenig 
beachtet und nur in einer Fußnote erwähnt. Erst ein vor einiger Zeit erfolgter freundlicher Hinweis von Herrn Prof. Dr: 
L. Ooliatz, Hannover, der die Formel 3. Ordnung in der genaueren Form mit Vorteil angewandt hatte, machte mich 
auf den praktischen Nutzen und damit zugleich auf den ganzen hier ängeschnittenen Fragenkreis aufmerksam. 


nung im Vorteil bleibt, sondern daß es sic 
rentialgleichungen noch höherer Ordnung zu werfolgen. 


Im folgenden wird die Erweiterung des Runge-Kutta- Verfahrens auf Differential- >: =” 


gleichungen beliebig hoher Ordnung durchgeführt, wobei freilich Gleichungen von höherer als 
der vierten Ordnung praktisch nur selten interessieren werden. Sowohl für das eigentliche 
Runge-Kutta- Schema, die „‚Schrittrechnung“, als auch für die Runge-Kuttaschen 
Mittelbildungen läßt sich dabei ein allgemeines-verhältnismäßig einfaches Bildungsgesetz angeben. 
Der Vorteil der so gewonnenen Formeln gegenüber der Behandlung der Differentialgleichung 
durch Aufspalten in ein System von Gleichungen niederer Ordnung liegt einerseits in der mit der 
Ordnungszahl n ansteigenden Genauigkeitsstufe, andererseits in der Reduzierung der Zeilenzahl 
des Runge-Kutta- Schrittes von vier auf drei in einem wichtigen Sonderfall. 


2. Rechnungsgang des allgemeinen Runge-Kutta-Verfahrens. 
Gegeben sei eine Differentialgleichung n-ter Ordnung - 


ya EP N N ee RE (1) 
mit den Anfangsbedingungen 
ya yon YEy Yay an MIET gerne (2). 
Gesucht werden Näherungswerte yj, Ya» = +5 Wr War - +3 Wr Yas see i neu; MTD, ul, 
an den Stellen 27 =&%,+h, 23 = x, 42h, ... nach einem Integrationsverfahren, welches als 
Verallgemeinerung des Runge-Kutta- Verfahrens angesehen werden kann. - 
Der Runge-Kutta-Rechnung legen wir im folgenden ein bestimmtes Schema zu- 
grunde, welches durch die gewöhnliche Runge-Kutta- Rechnung und insbesondere durch 
das oben erwähnte Nyström- Verfahren?) nahegelegt wird. Die Rechnung besteht dabei 
wie üblich aus drei Teilvorgängen, nämlich einer vorbereitenden ‚‚Schrittrechnung‘“, einer ‚„‚Mittel- 
bildung“ und schließlich der ‚‚Endrechnung‘‘ zur Bestimmung der gesuchten Näherungswerte 
Y» Ys -.., y"—D am Ende des ersten Runge-Kutta- Schrittes. 
Die vorbereitende Schrittrechnung vollzieht sich, wenn wir zunächst von 
einer Vereinfachung in Sonderfällen absehen, in vier Zeilen I, II, III, IV an den Argumentstellen 


h | 
a an DE Balae 1 =%77:; zy=®oth. 


Zu diesen Stellen gehören die vier Wertesätze y;, Yi, Yi, -- -, y" ® miti=[I, II, III, IV. Die 
Werte der ersten Zeile I sind die gegebenen Anfangswerte (2); die der drei übrigen Zeilen sind 
Hilfswerte. vorläufiger Art. Diese werden formal als Taylor- Entwicklungen von der Aus- 
gangsstelle x, aus und mit dem Argumentschritt h/2 für die Zeilen II und III, dem Schritt A für 
die Zeile IV gebildet. Sie benutzen demnach in ihren Gliedern die gegebenen Anfangswerte %,, 
Yo Yor » +» Ye" und fügen als letztes ein Glied mit y hinzu, welches gemäß y” =f(z, y, y', 
..., 49") mit Hilfe der Differentialgleichung berechnet wird, jedoch nicht mehr durchweg 
für die Ausgangswerte x, Y% Y,, - --, y®—D, sondern für die Werte an einer der Stellen I, II 
oder III, wie im einzelnen noch anzugeben ist. 

Wir erläutern den Aufbau der Taylor- Ausdrücke an einem Beispiel. Bezeichnen wir 
die für eine der Stellen © = I bis IV aus der Differentialgleichung berechneten n-ten Ableitung y® 
kurz mit f;, also: 

(a, yı, yh, Ye ya fi G=LFILUE IV) TEE aa or = Feat, (3), 
so nehmen in dem von uns gewählten Schema der Runge-Kuttaschen Schrittrechnung 
die oben erwähnten Hilfswerte etwa der Stelle II als Taylor- Ausdrücke mit dem Argument- 
schritt h/2 die folgende Form an: 


2 h n h? ZZ h3 hrz=1 hr 
un on hftpem 
1, ’ " h 22 h? oT hr -2 hr—1 
munter "mm gl ramı (n — 1)! 
h Ar—3 A— 2 . 4). 


’ „ 1777 ER 
Yı = % + 5) +...+y# N 


aeg 


Der See Be a U a TE 
4... 0 ee jet. a (a ie Den Se DER EEE Eee 


AR er u Beh 
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Für die Zahlenrechnung ist es zweckmäßiger, alle Ableitungen: y'® dimensionsgleich mit 


y zu schreiben, und zwar in der Form y®:h@/g!, in der sie in der ersten der Gln. (4) auftreten. Wir 
führen dafür die Abkürzungen 


ha 


Se (=1,2,...,n—1) OT ARE ET AR ER (5) 
ein und schreiben ähnlich für das letzte Taylorglied die Abkürzung 
hr - . 
ki lee (& Ya din don or, ) G=LIEUOLIV)....,. (6). 


In der Zahlenrechnung benutzen wir also die Differentialgleichung ( l1)in der abgeänderten Form 
(6), in der alle Ableitungen y;® in der Form v,; bzw. k; auftreten 5). Indem wir so ausschließlich 


mit y dimensionsgleiche Werte verwenden, erhalten die Taylor ausdrücke (4) folgende end- 
gültige Form: 


1 1 1 1 1 

Yır = %o +5 Yon Yo za vo-t. Sr are Un— 1,0 ron 

1 1 1 BR 1 
rn vo Bon (7) Yo ne (1) Yo ERS A % 1 ) Un— 1,0 + mi (1) S 

| e 1 mo1 1 (7). 
u tz ut tl te 

1 N 

%—10 = +5, ı) 5 


Für diese und weitere entsprechende Ta ylorausdrücke der Hilfswerte bedienen wir uns 
zweckmäßig einfacher Symbole. Wir schreiben für einen Taylorausdruck der Größe v, mit 


1 
dem Endglied k; bei halbem bzw. ganzem Argumentschritt }/2 bzw. h die Kürzung T, > %) bzw. 


T,(1, k;), wobei für g=(0 an Stelle von v, der Wert y tritt, also: 2 

2; > k) bzw. T, (1, %,) für die Taylorausdrücke von y , 
2 

T, nn k) bzw. T,(1,k) für die Taylorausdrücke von v,, 
2 


1 = ; k) bzw. 7, (1,.%k,) für die Taylorausdrücke von ®,, 
2 
DS k) bzw. Tu_1(1,%;) für die Taylorausdrücke von v,_ı, 


und es ist 


1 1 1 1 3 
15 =} )®% + E Jrer2o+ eh 
1 n—1 m 
:E Ing-i | q ) Un— 1,0 er (2) k; 


T,(1, k,) = +? ) Vg + lt) %+20..- 


ee he 


Mit diesen Kürzungen lautet nun die von uns zugrunde gelegte Vorschrift der Schritt- 
rechnung des allgemeinen Runge-Kutta- Verfahrens folgendermaßen: 


nen erget): 


5) Vgl. etwa die Zahlenbeispiele der oben angeführten. Arbeit 2), S.112 und 114. 
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Ber Fa Le a 


RER A ei ei; 
u 2 4 % MENT. —. 


Er; N 


sr Y-h ..- 


TEE Urzetıs 


Un — 2,1 U%n— 2,0 > “ wg 
A Me 41: 
n—11=%n—1,0 Se EN mn Tu— 3°) e 
Bi. k =k(a, 9 Yun ++ dm) Kr klaı Yın-Parpe -er> m; {E 
ae j . ? NT GE a © 
: 5 ‚Zeile I: 2 =, Zeile IV; av=ath 5 ir 
Be - Ynı = Yır Yıy = T,(l,kım) 
Bar j 2 h ns 
Br: N % m Pın Eu yy=T, (km) 
x | ine m 2,11 A %— 3, = Tn—2 (1, kın) 
a De : k h : 1 L nu 
2 | %—1,11= Tn—ı (5: kn) %”—1,w = Ta—ı (1, km) 
En: j' 2 . Q 
F; kıı = k (zır, Yu «+ 9%n—1,10)5 kıy = k (av, Yıys Brvs +++, 901,1) 
= R ‚Als Rechenschema nimmt das Runge-Kutta- Verfahren die folgende Form än: ri 
3 Tafel 1. Schema dos Runge-Kutta- Verfahrens für Ditferentialgleichungen n-ter Ordnung, 
3 i | x | Y | ra | Ye | ME oT | k; 


h I %o Yo v0 
i = 
u |s+4z|m, = $ %) T, Ex e) 


h ne ae 1 
„IV %+ h T, (, kın) T, (l, kın) 


Der Rechnungsgang ist somit: 
la) Eintragen der gegebenen Anfangswerte x Y9 %ıo ».-, ?ua—ı,0 in Zeile I; 
1b) Berechnen des zugehörigen Wertes k; nach der Differentialgleichung; 
2a) Bilden der Taylorausdrücke für Yy7 dırp »--» %a—ı,ır in Zeile II mit kr; 
2b) Berechnen des zugehörigen Wertes kıı nach der Differentialgleichung; 
3a) Wiederholen der Werte yır, %ıın =.., %n—2,ır in Zeile III mit Ausnahme des letzten 


N 
Wertes %— 1,1 = Tai WW kır = Mm—1,0 Bo kr; 
ge 


3b) Berechnen des zugehörigen Wertes kırı nach der Differentialgleichung; 
4a) Bilden der Taylorausdrücke für yıy, vılv; »+., %a—ı,ıv, in Zeile IV mit kur; 
4b) Berechnen des zugehörigen Wertes kıy nach der Differentialgleichung. nr 


Damit ist die vorläufige Schrittrechnung abgeschlossen. Esfolgt die Bildung vonnMittel- 
werten‘, k',k",..., k"— aus den vier Werten, ky, kır, ırr, kıy nach besonderen noch zu 
untersuchenden Vorschriften. Im Schema der Tafel 1 ist eine fünfte Zeile zur Aufnahme dieser 
Mittelwerte vorgesehen. Mit ihnen erfolgt dann die Berechnung der gesuchten Endwerte 
Y1 dın Pan +++, %n—ı,ı Zum Argument x, = x,—+-h, und zwar wieder in Form von Ta ylor- 
ausdrücken mit ganzem Argumentschritt, welche sich von den oben angegebenen Ausdrücken 
T,(1,%;) nur dadurch unterscheiden, daß das letzte Glied mit k; ersetzt wird durch den ent- 


» 
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sprechenden Mittelwert %', und zwar ohne weiteren Zahlenfaktor. Schreiben wir für den 
Taylorausdruck T, (1, k,) 0 hne das letzte Glied mit k; die Abkürzung T, (1), also 


ee, 


mit 990 = Yyofür q = 0; 50 lautet die Endrechnung: 


2 
)+R" 


Diese Größen bilden die Ausgangswerte des anschließenden neuen Runge-Kutta- Schrittes. 


3. Bestimmung der Mittelwerte kin, 


Aus den vier Werten kr, kır, kır, kıy sind nunmehr n Mittelwerte k, k’, k”’, ..., k"—D der- 
art zu bilden, daß die hiermit nach (11) berechneten Endwerte Yard, Vnidt Br den Werten 
Y1» %1 ++» ®n—ı,, einer Taylorentwicklung in möglichst vielen Taylor gliedern überein- 

. stimmen, ein Prinzip, welches bekanntlich auch zur Aufstellung des gewöhnlichen Runge- 
Kutta- Verfahrens geführt hat. Man geht somit ganz ähnlich wie dort vor und bildet: 


1. die Taylorentwicklungen der Werte y,=y(xz,-+Äh), %ı =tı (a, +h 
Zi Un— ni (&o — h) s > 
2. die Taylorentwicklungen der vier Werte %, kr, kım, kw: 


)» 50 n—11= 


Die erstgenannten Taylo r entwicklungen lauten: 


E hr min. hn+2 
Fe ren re 
+) eh) Hl wre ns 
=.) +) De )uetn rl wre + “3, 
a me) mir +21): en 
n +2 
He et 


Hierin enthalten die Ausdrücke 7, (1) ausschließlich die bekannten Anfangswerte (2). Im nächst 
folgenden Glied ist 4 aus diesen Anfangswerten mit Hilfe der Differentialgleichung (1) als 
'Funktionswert f berechenbar. Die folgenden Ableitungen yrtD, yr+2, ... sind dagegen aus 
durch ‚Weiterdifierenzieren zu bilden, und da feine Funktion der n Ka 1 unabhängigen Variablen 
2% YYıs 2.2. yD ist, so führt dieses Weiterdifferenzieren zu einer sehr rasch ansteigenden 
Anzahl von Einzelausdrücken, welche man zweckmäßig in gewissen Kürzungen zusammenfaßt, 
ähnlich wie dies bei der Herleitung des gewöhnlichen Runge-Kutta- Verfahrens üblich ist. 
Wir verwenden die folgenden Bezeichnungen: 


y®—u, =1,2,..0—1), et 


of BE (13) 


ya FRE 2,01), 9x — fh, 


sowie einen Differentialoperator D, der, angewandt auf eine beliebige Funktion 9 =9 (%, Y, U, 
.., %n— ı), die folgende Bedeutung hat: 


Do=9 Fu t%Ppı tr... Fun —at/pm—ı: nn ne (14). 


v r ER = & Be 
a el Er a IE Er 
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Ferner "bedeutet: : 


a RER er ER ) ee 
et tet Fe pl 3 | nt $ 


ne 
ee. 


Darüber Ne fassen wir noch einige ständig wiederkehrende Verbindungen zu weiteren Kür- 
zungen zusammen, nämlich: 


M=Dftwftwht.. the  * | | Is) 
Dif= Df4+3wDfo tw Dfit+.. HD )tuftuht. His), 


Hiermit lassen sich die höheren Ableitungen über y'” hinaus in bequemer Form anschreiben. 
Man erhält für die vier ersten Größen: 


y nk e 

y"tnD=Df u DALE) 
ya+9=Dif+fn—ıDf 

yr +9 —Dif+fn—ıDif+3 Dfa-ıDf+H(fa— tn) Pf 


Der Aufbau und die Gliederzahl dieser Ableitungen ist somit von der Ordnungszahl n der Diffe- 
rentialgleichung völlig unabhängig, ein Verhalten, welches sich auch bei den 1 aylorentwick- RT, 
lungen der k;-Werte wiederholt. Diesem Umstand ist es zu verdanken, daß sich das.Runge- 
Kutta- Verfahren ohne besondere Mühe auf Differentialgleichungen beliebig hoher Ordnung 
verallgemeinern läßt, ja daß sich für die gesuchten Mittelwerte sogar ein allgemeines und ver- - 
hältnismäßig einfaches Bildungsgesetz ergibt, wie wir noch sehen werden. 

Die Taylorentwicklung der vier Werte kı bis kıy führt auf die folgenden Ausdrücke, 
wobei wir künftig den Index O0 .zur Kennzeichnung der Ausgangsstelle fortlassen: 


E 
j 
“ 
{ 


h 
N 3 

BEE BIER CE a LER Er Do 
kur ae Te a 

Bei 1 mean An +2 
ST ea ee ira 18). 

1 hr +3 
ra! if +3f-ıDif+6DhmıDN——- + +: 

h" hr+l 1 hn+2 

kıy IA +h-ıDf) = 
LI®. 3 > Be mar hu+3 
+ [Pit DI +3 Da DI+ ZRH Dt: 


Aus den vier k,-Werten sind n Mittelwerte %® nach der Vorschrift 


ka) =Yıl kr Ya krı Yo kr -H-Yyıv kıv (q —ttE 1: PR Ey | Yunanı 1) ey (19) 


zu bilden derart, daß die mit ihnen gerechneten Endwerte (11) mit den Taylor werten (12) 
in möglichst vielen A-Potenzen übereinstimmen. Setzt man in (11) die Ableitungsausdrücke (17) 
ein, so erhält man die Forderung 


o—_ 2) x ee) HER % ” 3 RE 

vo (r)y2 +(”7 Diary JBirtnDn + @0: 
Ein Abgleich der beiden Ausdrücke (19) und (20) für beliebiges A fordert Gleichheit in jeder einzel- 
nen h-Potenz. Dies ist natürlich nur für endlich viele A-Potenzen erfüllbar, und da.den einzelnen 
h-Potenzen eine ständig wachsende Anzahl von Einzelgliedern zugeordnet ist, so ist dem an- 
gestrebten Abgleich von Mittelwert und Taylor wert bald eine Grenze gesetzt, eine schon 
vom gewöhnlichen Runge-Kutta- Verfahren her bekannte Erscheinung. 

Zur Bestimmung der vier Konstanten 
entsprechend den rechts in .Gleichung (20) 
den Abgleich der Glieder bis A" +2 
folgende vier linearen Gleichungen: 


Yaı Dis Yyrv stehen gerade vier Bedingungsgleichungen 
aufgeführten vier Gliedern zur Verfügung, wenn man 
verlangt. Hierbei ergeben sich, wie man leicht nachprüft, 


E32 e h 3 any { Es E= 
n! | trat. au Yav = (a) ne 


F 3 gran +4 5, Yan > | E22, = 93 I, 


245 a ; % 


gr en gran + 2 Dan | 
& = ee Je | =? En TR: a 
„a ER Lösungen erhält Ran nach einfacher rechnung er rechten. ‚Seiten ae P= = ng: 
N er RE Hirzaıı © ERROR EN 
Ar ü Be ; f . ; 3» 1 i E - = es 
Re ee Kr 2 armer En ® : (22% | 
ar RER Aresın |. An: 


Damit it sich für die gesuchten Mittelwerte die Re Formel: 


2 "= Gosmormter hand Hu + a 2 


Be 422. Dieser Kledtack der Mittehwerfe ‚enthält die geh n R Ditterentialgleichung Ba 
$ lediglich als Vorfaktor. Der eigentliche Aufbau der Mittelbildung hängt nur von der Zahl p ab, j 
By um welche:die Strichzahl q.des Wertes k(® gegenüber der Ordnungszahl n der Differentialgleichung es: 
erniedrigt ist. Damit weisen alle Mittelbildungen der höchsten, der zweithöchsten usw. Strich- 
zahl bis auf den erwähnten Vorfaktor für Differentialgleichungen aller Ordnungen n den gleichen 
Aufbau auf. Mit erhöhter Ordnung treten lediglich neue Mittelbildungen zu den im wesentlichen 
. unveränderten alten hinzu. In Tafel 2 sind die Zahlenwerte der Mittelbildungen nebst Vorfaktoren 
für n=1 bis 10 zusammengestellt, wobei freilich so hohe Ordnungszahlen praktisch nur selten 
interessieren werden. Der bekännte Mittelwert des gewöhnlichen Runge-Kutta-Ver- 
fahrens sowie die beiden Mittelwerte des N-yström- Verfahrens sind natürlich als Sonder- . 
fälle unter n—=1 und 2 mit enthalten. 


Tafel2. Aufbau der Mittelwerte k(”—»' in Abhängigkeit von der Ordnungszahl n der Differentialgleichung | 


und von p für n,p=]1 bis 10. 
{ Feste Faktoren von Veränderliche Vorfaktoren für n = 
a ISBN ZERTEEt 
Sr sEreen 2 
1 1 43 1 1/6 |1/73 |ı/2 |2/3 | 5/6 1.776 © 478-|-s3ja9lsa > 3 
2 Az 4 0 vı2 |ıJ7a [172 1876 |.5/4 | 7a | 78.» 3 15/4 : 
3 9 BR > 1720 11/5 |172 ı |77a |145 21/5 | 6 2 
4 16 8 2 130 1.176: 1 172° 9/6. „u 22]3. 91/8 v7 a 
< 5 25 10 —3 1/42 | 1/7 1/2 4/3 3 6 ER 
6 36 12291 1/56 | 1/8.:| 1/2 | 3/2 |15/4 
ae Ag Er Tl er 1/72. | 1/9 | 1/2 | 5/3 2 
= 64 16 —6 .1/90| 1/10 | 1/2 Mn 
9 81 18 ES; 1/7110 1/11 : 
10 100 20 8 ß 1/132 


4. Die Bas averbiltulse. 

- Mit den soeben festgelegten Mittelwerten k'® wird, wie aus ihrer Herleitung folgt, ein Ab- 
gleich aller Größen v, = y'® ha/q! bis zur Potenz h" +? einschließlich erreicht. Die bedeutet 
eine Genauigkeit der Ordnung h" +2 für y, h"+1für y’ usw. bis A? für y"7V als letzte der zu 
berechnenden Unbekannten. Bei einer Behandlung der Differentialgleichung nach dem gewöhn- 
lichen Runge-Kutta- Verfahren durch Aufspalten in ein System von Gleichungen. 1. Ord- 
nung erreichen demgegenüber sämtliche Unbekannten die Genauigkeitsordnung Ah*. Soll also das 
unmittelbare Integrationsverfahren mit dem alten Runge-Kutta- Verfahren wettbewerbs- 
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fähig sein, so wäre eine Genauigkeitsordnung von h* auch für die Größe yır=» erforderlich. Hier- 
zu ist aber für yr—D, also für 9, ı ein Fehlerabgleich auch noch in den Gliedern der Ableitung 
y@ +3) notwendig, was einen Abgleich in dem folgenden vier zusätzlichen Gliedern 'E 


Dif, fn—ıDif, 3Dfn-ıDf und (fi-ıtfn—)Df 


der Ausdrücke (17) und (18) voraussetzt. ’ 
Die Nachrechnung zeigt nun, daß dieser Abgleich durch die aus den früheren schwächeren 


Forderungen bestimmten Koeffizienten = 


Y PLE- 


n n . ER 
YAaA—- LIT -ıW Tr Yn— 1,1 = Yn— 1,11 7 } ‚ E: i 
von selbst mit erfüllt wird. "Durch diesen Umstand, der übrigens durch die hier im wesentlichen 


vom alten Runge-Kutta- Verfahren her übernommene Art der Schrittrechnung herbei- 
geführt wird, ist die Mindestgenauigkeit der Ordnung A? für alle zu berechnenden Größen ge- 
sichert. Es ist somit 


y'® genau von der Ordnung h"=@ +? für n—q 22, 
y“) genau von der Ordnung A für n—q=1. | - 


In Tafel 3 haben wir die Genauigkeitsordnungen noch einmal für n=1 bis 5 im einzelnen an- 
gegeben, wobei n=1 dem alten Runge-Kutta-Verfahren, n=2 dem Nyström- 
Verfahren entspricht. Es zeigt sich, daß von n=3 an das hier beschriebene unmittelbare 
Runge-Kutta- Verfahren für Differentialgleichungen höherer Ordnung einem Aufspalten 
der Gleichungen und Behandeln nach dm Runge-Kutta- Verfahren 1, Ordnung oder 
dem Nyström- Verfahren hinsichtlich des Genauigkeitsverhaltens eindeutig überlegen ist. 


Tafel 3. Genauigkeitsordnungen des unmittelbaren 
Runge-Kutta- Verfahrens für Differential- 
gleichungen n-ter Ordnung bei n=1 bis 5. 


q | y | RE 1] 2 | 3 | 4 | 5 
0 Yy h% ha h5 ns a? 
1 y pr Top Iope re 
2 gi | | 3 
3 772 hi h* 
4 el hi 


5. Arbeitsersparnis in Sonderfällen. 
Bei Differentialgleichungen n-ter Ordnung hat mian es häufig nicht mit dem allgemeinen 
Fall einer Abhängigkeit der Größe y'* von sämtlichen niederen Ableitungen y® mit g=0 bis 
n—1 zu tun. Fehlen nun in der Funktion f eine oder mehrere der Größen y®, so brauchen die 
betreffenden Hilfsgrößen yP, yi}, y\9) bzw. vyrr, vor, vv in der Schrittrechnung nicht gebildet 
zu werden, während der Hauptwert y!2 bzw. v,ı der ersten Zeile als Ergebnisgröße natürlich im 
allgemeinen wohl mitgeführt wird. 


Fehlt insbesondere in der Differentialgleichung eine Abhängigkeit von der zweithöchsten 
Ableitung „=D (was ja z.B. bei linearen Differentialgleichungen durch eine bekannte Um- 
formung stets herbeigeführt werden kann), so fallen die beiden Zeilen II und III in eine einzige 
Zeile II/III zusammen, und die Schrittrechnung reduziert sich von vier auf drei Zeilen, wobei 
in den Mittelbildungen der mittlere Wert kyıyıır jetzt doppelt zu werten ist. Das bedeutet gegen- 
über dem alten Runge-Kutta- Verfahren eine mitunter wesentliche Arbeitsersparnis, da 
in der Berechnung der Funktionswerte %; nicht selten der weitaus größte Arbeitsanteil steckt; 
vgl. hierzu auch die Zahlenbeispiele in ?) sowie das Beispiel am Schluß dieses Aufsatzes. 


6. Das Verfahren für Differentialgleichungen 3. und 4. Ordnung. 

Zum Schluß sei das in Tafel 1 für den Fall der Differentialgleichung n-ter Ordnung in all- 
gemeiner Form angeführte Rechenschema für die praktisch wichtigsten Fälle der Differential- 
gleichungen 3. und 4. Ordnung noch einmal im einzelnen wiedergegeben; bezüglich des Ny- 
ström- Verfahrens für die noch wichtigere Gleichung 2. Ordnung vgl. 2). Indem man die Tay- 


= 1 
l or ausdrücke 7, 5 k) und 7, (1, %,) der Tafel 1 nach Formel (8) für n = 3 und 4 ausschreibt, 
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erhält man die beiden Schemata der Tafeln 4 und 5, wobei für &, jetzt einfachheitshalber 0 ge- 
setzt ist. Die Mittelbildungen sind für n=3 und 4 der Tafel zu entnehmen; hier ist, soweit 
erforderlich, noch eine für die Maschinenrechnung PaIanTE a hinzugefügt worden. 


were Ordnung: 
"= y,y") 
’ m’® 
5 =" hm yonv,) 


Tafel 4. Schema.der Runge-Kutta- Rechnung für Differentialgleichungen 3. Er 


» ‚h? 

s | ® | y v=yh a k; 
DU | Yo u2 vo v0 kı 5: 

1 

u h/2 er Del rn 4 ?20 SZ sh vYot ot = Ch dot > SE kır 2 
II h/2 ” ” ” „ ” Vo + 3 = n k;ıı 
IV h Y%H+ vo + do + kıt vt2vYo + 3 ku v%0+ a kıv 
Mittel: k | 2 ee k” | 

I | h | Yı=YtVtotVtaortk | a U men Re Y= to + 

Mittelbildung: 
k ni Ikr-+6kır +6 kın — kıv) = 0,45 kı 0,3 kıı —0,3 kıır — 0,05 kıv, 
—=kı +kı + kın, 


I == a (kı +2 kır +2 ka + kıv) == 0,5 kı —+kıı —+ kyır 0,5 kıv. 


Differentialgleichung 4. Ordnung: 
yY=f(ayy\y,y”) 


k; Ben %,9, 01,03, 05) 
ı Y 54 (%, % %1> Ua, da 


Tafel 5 Schema der Runge-Kutta- Rechnung für Differentialgleichungen 4. Ordnung. 


[73 h? [773 h3 
i|® Yy v=y'h »=Yy > 3=Y u 
I|0 Yo do Vo | Vzo | kı 
| j 
| 1 1 3 3, | 
11 = Yo+ 3ri0+ 4 ro tg tot jgkt vortoary En hy %0 +5 90 +ty%k Yo t+2k; | kı 
] 
III z ” „ ER) ”» ’ 27 2 > EZ) ort 2k,ı kın 
IV Ah yo+ do + %o + Yo + krıı | Yo +29 + 390 + Akıyı \®go + 3.090 + 6kım) Yo + #kın | kıv 
Mittel: 2 k k’ k' Be, 


I | h | Yı = Yo + vıo + Paot do + k Win =Vıo+ 2% + 39904 Kdaı = dot 39a0+ kw, —=tgo+k”” 


Mittelbildung: 
1 
k - (8 kı +4 kıı +4 kın — kıy) = 5 (1,6 kr +0,8 kır +0,8 kın — 0,2 kıv), 
k' =+ (Ik +6 kı +6 km —kıy) =1,8Kı —1,2 kır + 1,2 kın — 0,2 kıv, 


= = +2kı +2 kıı 
z 
en (kr +2kı +2 km Fkyv)=— (2kı +4kı +4 kım +2kıy). 
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7 Zahlenheispiel einer Ditferentialgleiehung 4. Ordnung. SR 

- Das Verfähren für Differentialgleichungen 4. Ordnung sei an einem € Zahlenbeispiel 
erläutert. Gegeben sei die Differentialgleichung _ | mE ER 
. 2 - a ne 

mit den Anfangsbedingungn y=y u = I für ,—0, wobei die exakte Lösung _ ee 
mity=e* en ist. Die Schrittrechnung reduziert sich hier auf nur 3 Zeilen, und die Zwischen- Se 
werte für v,, v, und v, brauchen nicht mitgeführt zu werden, Tafel 6. Die Stellenzahl in den. 


Tatfel6. Runge-Kutta- Rechnung der Differentialgleichung y!Y = y. Schrittgröße h = 0,2. 


i | © | er % u v3 ki 
1,0000. 0000 | 0,2000.00000 | 0,0200 0000.00 | 0,0013 3333 333| 0,0000 6666 6667 
IL/IIL 0,1 1,1051 7083 -7367 8055 4 
iv 0,2 1,2214 0701 8142 7134 
Mittel 0,0000 69422 0,0002 805419 
I 02 | 1,22140276 |’0,244280542 | 0,0244 2804 56 | 0,0016 2853 734| 0,0000 8142 6851 
I1/LIT 0,3. „| 1,3498 5870 - 8999 0580 
IV 0,4 1,4918 2988 9945 5325 
- Mittel 0,0000 84792 0,0003 426547 | 0,0005 2281 602| 0,0003 6056 300 
0,4 1,4918 2468 0,2983 64900 0,0298 3647 37 | 0,0019 8910 034 


einzelnen Spalten wird so gewählt, daß überall etwa die gleiche Anzahl geltender Stellen 
auftritt. Die dadurch von Spalte zu Spalte um eins vermehrte Dezimalenzahl empfiehlt sich 
auch mit Rücksicht auf Abründungsfehler, welche sich sonst durch die Faktoren der Schrittrech- 
nung und Mittelbildung leicht unzulässig vergrößern können, 


Tafel7. Fehler der Größen %, y’, y”’, y'’’ in Einheiten der 6. Dezimale für das Beispiel yIY = y; = ef. 
RK=Runge-Kutta. 


RK-Verfahren für DGn. 4. Ordnung Nyström-Verfahren | pp.weeehren 
% — un 

’ 073 1273 ’ ’ [273 Y == y’ er 
Yan Yy Y Yy =y" = [7 [7 

| | | us v-y vr yıı 

0,2 0,00 —0,05 —0,48 0,25 — 2,76 0,56 — 2,76 
0.28 177°0,08 0,20 —1,02 0,56 — 6,07 0,74 — 6,74 
0,6 —0,08 —0,46 — 1,62 0.94 —10,15 0,39 — 11,34 
0,8 0,21 —0,85 —2,31 1,37 —15,30 | —0,60 — 20,10 
1,0 —0,43 —1,40 1 1,87 2,87 | 8 — 30,69 


In Tafel 7 sind die Fehler der berechneten Größen y, %', y', y'’’ für die fünf ersten 
Runge-Kutta- Schritte angegeben. Zum Vergleich sind die entsprechenden Fehler einer 
Rechnung nach dem Nyström- Verfahren und einer nach dem gewöhnlichen Runge- 
Kutta- Verfahren aufgeführt, wobei die Differentialgleichung in das System y’=u, u’ = y 
bzw. Y=u, wW=», v’ =w, w' = y aufgespalten ist. Man erkennt den großen Genauigkeits- 
vorsprung des neuen Verfahrens. 


Eingegangen 10, Dez. 1947. 
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Fehlertrennung durch Wahrscheinlichkeitsbetrachtung. 


Er aM be ung der Universitätssternwarte Bonn.) 
ER / Von Joseph Meurer. rs 


‚Es wird eine Methode angegeben, die Fehler zweier voneinander verschiedener Meßweisen ein und. der- 


PEST selben ‚Größe zu irenmen, wenn keine Anhaltspunkte dafür vorliegen, wie sich in der Praxis der beiden Meß- 

NE arten die Fehler in den Meßresultaten auswirken. Durch eine Analyse der Wahrscheinlichkeiten in der Ver- 
0 teilung der Differenzen, welche die Messungen der betreffenden Größe nach beiden Methoden ergeben, über 

== ‚die zu messenden Argumente der beiden Meßweisen ist es möglich, diejenigen Wertintervalle für die zu messende 


° BE zu finden, in denen sich die eine bzw. die andere Meßmethode besonders durch ihre Fehler bemerkbar 
2 A method is given to separate Ihe errors of two different ways of measuring the same quantity, if the 
‚effect of the errors on the results when practising the two ways of measurement, is unknown. The author 

analyzes the probabilities appearing in the distribution of the differences obtained by measuring the quantity 


particular influence. 


B »Toä pabore gaeTcd MeToN BbINeNIeHUN ONIMÖOR ABYX P&SIMUHBIX CIIOCO6OB H3MEpeHHun 
'ORHOH U TOÜ-ke BEJIHYHHBI U PASHOCTH Meskıy pesy.ıpTaTamn 000nx usMepennä ua cayuas, 
EeCHH AA 9TUX OÖOUX CIOCODOB HEHSBECTHO BINAHNE OMMÖOR HA PesyIsTar usMepennz. 
Pacnpegesaa pasHocru Meray pesyısTarTama 000Nx HSMepeHnf B 3ABUCHMOCTU OT W3Me- 


pAeMbIX IIO O00HM CIOCODaM &APIyMeHTOB U AHANNSUPyA BEPOATHOCTE ITUX pasHocrel,” 


MOKHO ONPeNelHTb Te HHTEPBAJEL APTYMEHTOB, B KOTOPEIX TOT HAIH Apyroi MeTon naer 
. 0CODEHHO\ DOMBILYIO OIMHORY. RE Er E 
DNA : 


I. 


Eine Größe H sei gemessen auf zwei völlig verschiedene und voneinander unabhängige 
Weisen. Die Bestimmungsstücke der ersten Methode, die zur Bestimmung von H gemessen werden 


. müssen, seien db, und 5,, diejenigen der zweiten Methode 5, ünd b,. Es werden dann nicht bei 


jeder Messung die dem tatsächlichen Wert H zukommenden b; (©= 1,2, 3,4) gemessen, sondern 
diese sind mit gewissen systematischen und zufälligen Fehlern ö 5; behaftet, welche an Stelle des 
wirklichen Wertes HZ bei jeder Methode einen von H verschiedenen Wert H, und H, ergeben, 
wobei im allgememen H, #H,. Es ergibt sich also für jede Messung ein Fehler A,, der zu H, 
gehört und einer A, für H,, die abhängen von b,, b,, öb,, öb, bzw. b,, b,, öb,, öb,. Es ist also 


H,=H-+4,(b,,b,, öb,,öb,) 1 
H,=H-+4,(b,,b, öb;, öb,) th: 


Es sei nun weiter angenommen, daß sehr viele gleichzeitige Messungen vorliegen für H nach beiden 
Methoden, so daß man also sehr viele zusammengehörige Wertepaare (H,, H,) hat, die sich alle 
auf voneinander verschiedene Werte von H beziehen sollen, wo aber jedes Paar zu ein und dem- 
selben Wert von H gehört. Man kann dann unter gewissen Voraussetzungen durch eine Wahr- 
scheinlichkeitsbetrachtung feststellen, wie sich praktisch die Fehler ö5, auf H auswirken. (Theo- 
retisch ist ja ihr Einfluß durch die Fehlertheorie der beiden Methoden gegeben, und die A; sind in 
ihrer funktionalen Abhängigkeit von b; und öb; bekannt.) Es ist dies vor allem dann vorteilhaft, 
wenn man nicht weiß, in welchem Umfang und unter welchen Bedingungen die ö5; bei einer oder 
bei beiden Methoden auftreten. An Stelle des Wertes H bekommt man die beiden Werte H,;, und 
es ist in einem solchen Falle zunächst nicht zu sagen, welche Methode nun den größeren Fehler 
enthält. 
Man kann aus allen vorliegenden, nach den beiden Methoden gefundenen Wertepaaren (4,, 
H,) eine Differenz D 
D=#,—-B,=|+#4,%.b)+4s(b,b)| . . .» 22... (2) 


: bilden, — die Abhängigkeit der A; von den öb; ist in (2) und im folgenden nicht mehr besonders _ 


angegeben —; für jede gleichzeitige Messung von H nach der einen und der anderen Methode 
erhält man dann einen Wert D. Gemäß (2) hängen die Werte D nur von A,und A, ab, den Fehlern 
der beiden voneinander unabhängigen und verschiedenen Meßmethoden. Es ist nun die Aufgabe, 
aus den Differenzen D die Fehler A; für sich auszusondern. 

In (2) hängen die Werte Du. a. von den d; ab, und diese d; sind als Bestimmungsstücke der 
Messungen bis auf die Fehler ö b; bekannt. Die 5; bestimmen maßgebend den theoretischen Fehler- 
verlauf, von ihnen hängt es ab, in welchen Intervallen von 5; die zugehörigen öb; besonders wirk- 
sam sind, d.h. die A; besonders große Werte annehmen, bzw. das Umgekehrte der Fallist. Hin- 
sichtlich D wirken bei der Abhängigkeit gemäß (2) von den 5; A, und A, zusammen. Es sei nun 
vorausgesetzt, daß die Messungen der Größe H so zahlreich sind, daß b, von b, und b, von b, als 
unabhängig angesehen werden können, daß also im Durchschnitt alle Kombinationen b,, b, bzw. 


 concerned by both methods, among the arguments of the two ways of measurement. T’hus, it is possible to 
find intervals of the measured quantity in which the errors of ihe one or the other method of measuring are of 


ar wi 
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b,, b, durchweg gleichhäufig sind und damit ebenso alle überhaupt nur möglichen, — durch die 


praktisch möglichen öb; bedingten —, Fehler A, und A, durchweg gleich oft vorkommen. Gemäß 


igkeit von ei indem man ein Streu- 
(2) kann man nun D in Abhängigkeit von einem der b,, etwa b,, darstellen, an 
ae zeichnet, dessen Abseisse’die gemessenen b, und dessen Ordinate die zugehörigen D-Werte 


enthält (Bild 1). In diesem Streufeld steckt dann der Verlauf von D mit A,, natürlich gestört 


Bild 1. 


durch die statistischen Schwankungen von A,. Nun sind A, und A, verschiedene Funktionen 
von b,, b, bzw. b,, b,. Je größer diese Verschiedenheit ist, je charakteristischer der Gang eines 4; 
mit einem 5; ist, um so mehr muß sich dies in einem solchen Streufeld bemerkbar machen, um so 
weniger wird er durch die Werte des anderen A; gestört. Wenn nun wegen der vorausgesetzten 
Unabhängigkeit der 5, auch die A; voneinander unabhängig sind und durchweg gleichhäufig vor- 
kommen, dann sollte es möglich sein, durch ein Ordnen der D-Werte nach einem.b; mit besonderer 
charakteristischer Abhängigkeit, was ja ein solches Streufeld ist, diejenigen Werte von D zu finden, 
bei denen das eine der A, besonders wirksam ist und bei welchen dies für das andere A; zutrifft. 
Es seien nun sehr viele vorliegende D-Werte nach 5, geordnet. Es sei 
N ß,),. = Anzahl aller D-Werte im Intervall $,<d,<Pß, 


N ö,,, = Anzahl aller. D-Werte im Intervall d,<d,<ß, mit DD, 


die Wahrscheinlichkeit dafür, daß ein beliebig herausgegriffener D-Wert, dessen 5, in. dem Inter- 
vall (ß,, ß3) liegt, >Ööist. P hängt gemäß (2) von Dab und natürlich auch von 6. Wenn nun 
jetzt A, eine sehr charakteristische Funktion von 5, ist und für A, obige Voraussetzung erfüllt 
ist, dann hängen diejenigen Werte von D, welche wesentlich durch A, beeinflußt sind, in ihrem 
Gang mit 5, ebenfalls von A, ab, d.h. sie weisen einen diesem ähnlichen Gang auf. Das muß dann 
aber auch für die Wahrscheinlichkeit P gelten. In einem Intervall (,, ß,), in dem A, groß ist, 
wird auch die Wahrscheinlichkeit groß sein, daß ein herausgegriffener D-Wert dieses Intervalls 
über einem vorgesehenen Wert ö liegt, jedenfalls wird sie größer sein als in einem anderen Inter- 
vall ßJ, <bd,<P, wo A, kleiner ist gemäß der theoretischen. Abhängigkeit von b,. Dabei ist 
natürlich wesentlich, daß der Gang des A, mit b, gemäß obiger Voraussetzung durch die statisti- 
schen Schwankungen von 4A, nicht unterdrückt wird. Man kann nun den Wahrscheinlichkeits- 
verlauf P mit b, für verschiedene ö, d.h. für verschiedene Intervalle der D-Werte, zeichnen in 
Abhängigkeit von b,, indem man für 5, hinreichend enge Intervalle (8, ß,) nimmt und P bestimmt. 
Für diejenigen Intervalle der D-Werte, die sehr stark von A, 
abhängen, sollte sich dann ein ähnlicher Verlauf von P mit 5, 
ergeben wie derjenige von A, mit d,. Dieser Wahrscheinlichkeitsverlauf kann 
dann in das Streufeld der D—b,-Wertepaare eingetragen und daraufhin geprüft werden, wieweit 
er in dem Streufeld ausgeprägt ist. 
IE 


Das Vorstehende soll nun an einem speziellen Beispiel erläutert werden. Es handelt sich 
um die Messung einer Höhe H über B (Bild 2) einmal durch Bestimmung der direkten Entfernung 


a 


Instrumentes immer wieder verändert wurde. (Diese Abhängigkeit ist im allgemeinen so, daß 


ge 3” ‚eine sehr ausgeprägte Fehlerkurve, welche sich wegen H =e-sin a darstellt als 


_ B, deren Methodik hier im einzelnen nicht interessiert, für die aber wesentlich war, daß die zu- 


Glied hohe Werte annimmt. Eine theoretische Fehlerkurve, die AH in Abhängigkeit von a dar- 


GE TR ME RER DB a ? et 
des Höhenwinkels a von A ‚aus, und zum anderen durch eine Messung von H direkt über 


fälligen und unsystematischen Fehler, — es handelte sich um die Messung mit Druckdosen —, . 


sehr stark überwogen, so daß notwendig die theoretische Fehlerabhängigkeit von H durch jene 
Fehler vollkommen überdeckt werden mußte, bzw. durch unkontrollierbare Beeinflussungen des 


_ mit zunehmendem H auch die Fehler A H größer werden.) Demgegenüber hat die erste Methode 


AH=4,=de-sina-+ecosada . . en EEE) 


. 3 E 


Ba Sie Bild 2a. 


In (4) ist jetzt in der Bezeichnungsweise von (1,574, 0b,= 00, bs = e und öb,=de. Aus 
(4) folgt nun, daß bei kleinen & das zweite Glied in (4) groß ist, während bei großen a das erste 


stellt, hätte also als Superposition der beiden Glieder in (4) einen Verlauf, den qualitativ Bild 2a 
wiedergibt, wobei der linke Teil der Kurve von dem 2. Glied und der rechte von dem 1. in (4) her- 
rührt. PURE : 
Der Fehler A, der zweiten Bestimmung von H direkt über B hat jene Winkelabhängigkeit 
nicht. Er kommt gleichmäßig bei allen & vor, ist über diese statistisch zufällig verteilt. _ 
Die D-Werte enthalten also gemäß (2) ein Glied, das sehr stark &-abhängig ist, nämlich A, 
gemäß (4), und ein zweites A,, das dies nicht ist. Darum können hier durch eine Wahrscheinlich- 
keitsbetrachtung von der vorhin skizzierten Art, — die dort gemachte Voraussetzung von der 
Unabhängigkeit der Bestimmungsstücke trifft auch hier zu —, die Fehler A, und A, in ihrem Ein- 
fluß auf H getrennt werden. 
Das bereits erwähnte Bild 1, aus technischen Gründen stark zusammengedrängt, und nur 
ein Ausschnitt des vollständigen Streufeldes darstellend, enthält D-Werte gemäß (2), wie sie sich 
aus den zusammengehörigen Werten HZ, und H, ergeben; es handelt sich um 560 Messungen, die 
in Abhängigkeit von dem jeweiligen gemessenen & dargestellt sind. Auf den Abszissen sind nicht 
die unmittelbaren D-Werte abgetragen, sondern ihre Prozente von dem gemessenen Wert H,. 
Das Streufeld wurde nun gemäß Formel (3) ausgezählt und P bestimmt für Intervalle von 5° 
zu 5°. Für ö wurde genommen 1%, 3% und 5%. Dem entsprechen die drei Kurven in Bild 3. 
Die oberste Kurve in Bild 3 stellt die Wahrscheinlichkeit P dar in Abhängigkeit von a, daß ein _ 
Punkt aus einem Intervall & beliebig herausgegriffen in Bild 1 über 1% liegt, oder daß ein D- 
‚ Wert des betreffenden Intervalls >1% des gemessenen H,-Wertes ist. Für D>3% gilt die mitt- r 
lere Kurve, für D>5% die unterste Kurve. Sie sind einmal ausgeglichen. Für die in Bild 3 an- ; 
„gegebenen «-Intervalle gilt genau immer der Punkt der Kurven, welcher in der Intervallmitte liegt. \ 
(So gibt z.B. die oberste Kurve an, daß die Wahrscheinlichkeit dafür, daß ein D-Wert im Intervall 
19°—24° > 1% ist, den Wert 0,77 hat.) Die Gestalt der Kurven ist nun wesentlich verschieden. 
. Während die Kurve für 1% nur schwach den Fehlerverlauf der Kurve Bild 2a aufweist, ist dies 
bei der für 3%, in sehr starkem Maße der Fall, während er bei 5% kaum noch ausgeprägt ist. Dies 
bedeutet jetzt, daß die Fehler der auf & und e basierenden Messung hauptsächlich 3% des H;- 
Wertes (d.h. in erster Näherung auch des wahren H-Wertes) ausmachen, während für die größeren 
D-Werte A, eine Rolle spielt, das die «-Abhängigkeit nicht hat, d.h. die Fehler A, überwiegen 
- die A, bei den größeren Abweichungen von H. Selbstverständlich gilt dies alles nur unter der 
vorhin gemachten Voraussetzung der Unabhängigkeit aller Bestimmungsstücke, d.h. ihrer tat- 
sächlich gefundenen Meßwerte voneinander, was bei den hier benutzten Werten im statistischen 
Durchschnitt zutrifft. 
Man kann die Methode noch weiter verwenden. Man kann annehmen, daß in jeder MeBß- 
serie die Fehler der Bestimmungsstücke innerhalb bestimmter Grenzen liegen, daß also in diesem 
Falle etwa |d« |< 1/10° und |de |<100 m ist. Ist nun ein durch & und e gemessener Wert HZ, 


vorgegeben, so muß der wahre Wert H in dem Intervall H,—4,<H<H,—-A, liegen, wo- 
bei A, durch die obere Fehlergrenze von 6a und ö e bestimmt ist. Nimmt man den zu 4, gehörigen 
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186 ’ Meurers, Fehlertrennung durch Wahrscheinlichkeitsbetrachtung er Ba 


Wert H,, so kann man D-Werte auf folgende Weise bilden: Die Werte sollen =D geseizt sein, ; 
wenn H,in das Intervall H, + 4, fällt. Wenn es außerhalb liegt, soll unter D immer seine Diffe- 
renz gegen H, — 4, verstanden sein, wenn H, < H,, und es soll diejenige gegen H, +4, unter. 
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Bild 3. "Bild 4. 


D verstanden sein, wenn H, > H,. Diese D-Werte, die man wieder in % von H, ausdrücken kann 
müßten durchschnittlich von den Fehlern A; frei sein, wenn die numerische Annahme für die 
obere Fehlergrenze von d& und de zutrifft. In Bild 4 sind, genau wie in Bild 3, die Wahr- 
scheinlichkeitskurven für die auf die beschriebene Weise erhaltenen D-Werte dargestellt, d.h. sie 
geben an, mit welcher Wahrscheinlichkeit ein beliebig herausgegriffener Fehler von diesen nach 
Abzug der Fehlergrenze von + 1/10° und + 100 m noch Übriggebliebenen >1%, >3% und >5% 
liegt (obere, mittlere, untere Kurve). Man sieht nun aus Bild 4, daß von 3% an der charak* 
teristische‘ Fehlerverlauf nicht mehr ausgeprägt ist, wie dies gemäß (2), (3) und (4) bei Bild 3 der 
Fallist; die nach Abzug der Fehlergrenze noch übrig bleibenden Fehler sind auf die kleinen %-te 
beschränkt; das bedeutet, die tatsächlichen Fehler sind noch etwas größer als obiger Fehler- 
grenze entspricht. 


04 


02 


Bild 5, 


Man kann die D-Werte auch in Abhängigkeit von e darstellen und erhält ein Streufeld analog 
dem in Bild 1. Für den Verlauf von D mit b, = e gemäß (2) und mithin auch (Bild 5) für den 
von P gemäß (3) ist wesentlich, daß vor jeder Messung eine Kalibrierung der Entfernyngsskalen 
vorgenommen wurde mit einem Testobjekt, das 14,0 km entfernt war.‘ Es ist also zu erwarten, 
daß die tatsächlich gemachten Fehler für diese Testentfernung besonders klein sind. In der Tat 
zeigen die Wahrscheinlichkeitskurven in Bild 5, die das gleiche bedeuten wie die in Bild 3 und 4 
für diese Entfernung mit geringen Abweichungen ein Minimum. Darin dürfte sich wohl die Emp- 
findlichkeit der oben beschriebenen Methode ausdrücken. 

Eingegangen 29. Dezember 1947. 


'erei VerfahrenzurSpannungsberech- 
nung. nnwandigen pri 
 körpern unter Innendruck. 


„In dünnwandigen prismatischen 
© sich bei Innendruck einerseits reine Längsspannungen 
0x aus (Bild 1), welche im allgemeinen als konstant an- 
2. genommen werden können undsich aus der Gleichung 


> errechnen. Hierin ist p der innere Überdruck, U der - 


\ 6 Bu | 
RETTEN { 


a 
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"Bild 1. Querschnitt einer dünnwandigen Röhre unter 
j Innendruck.. B 

F die von der Wandungsmittellinie eingeschlossene 
Fläche und h dieWandstärke. Andererseits treten Um- 
fangs- oder Ringspannungen o; auf, die von der Längs- 
kraft S und dem Biegemoment M (Bild 2) herrühren. 
Die Gleichgewichtsbetrachtung an einem durch zwei 
benachbarte Achsialschnitte begrenztern Wandungs-. 


streifen mit den rechtwinkligen Koordinaten yund 2 


führt auf drei Gleichungen, welche die Berechnung der 
. drei Spannungsresultanten (M Biegemoment, 8 Längs- 


f 


Bild 2. Kräftespiel an einem Wandungsstreifen.‘ 


kraft, Q@ Querkraft) des an beliebiger Stelle gedachten 
linken Begrenzungsschnittes erlauben, wenn die Größen 

" Mo Yo» Zo, bekannt sind; diese stellen die Spannungs- 
resultanten des als fest vorzustellenden rechten Begren- 
zungsschnittes dar (Koordinaten %,, 2), Biegemoment 
M,, Schnittkraft in y-Richtung Y,, Schnittkraft in z- 
Richtung Z,). Es ergibt sich 


E MITTEILUNGEN 
smatischen ‚Hohl- | 


Hohlkörpern bilden = 


Br ER Umfang (Abwicklungslänge der Wandungsmittelline), 


hr 


ersetzt: . VERASEEA ra - j 
= Zu R——l.. 
& Hn=P (0 — 2) $ ee: ERRER (8), 
Z=-Ph—Y). EEG Ba 2 Fahr (6). 
. : Nach Einsetzen in das Gleichungssystem (1), (2), (3) 
folgen: RR ; EN GER 


MM + Ne) HAN 


+ He] } TO 
S=[Yo+ pP @—a)] 008 a—[Z+P (y—ylsina (2), 
Q=[N+P @—z)]sina+lZ+P(y—yleosa (3). 
‚Die Aufgabe ist: daher im allgemeinen dreifach sta- 
tisch unbestimmt. Die drei statisch Unbestimmten 


M,, Y, und Z, werden zweckmäßig mittels der folgen- 
den Substitutionen durch die Längen y,, 2, und 


MP ywyP tere)... 
S=rle—z)esa + (y—y)snal..(&, 

Q = pl@— 2) sina — (y— y)cosal. . (. =. 

Aus dieser Darstellung geht unmittelbar folgender geo- 


metrischer Sachverhalt hervor: - r 
Es existiert ein besonderer Querschnittspunkt O,, 
der als „„Biegezentrum‘“ der dünnwandigen Röhre be- 
zeichnet sei; er besitztfolgende Eigenschaften (s.Bild 3): 
1. Das Quadrat des Abstandes vom Biegezentrum, 
vermindert um die Konstante 0°, ergibt nach Multi- 
plikation mit dem halben Betrage des Innendruckes 


Bild 3. Bestimmung: von Längs- und Querkraft mit Hilfe 
des Biegezentrums. 


das Biegemoment. Daher sind die Schnittpunkte 
des Kreises vom Radius o um das Biegezentrum 
mit der Wandungsmittellinie zugleich die Momen- 
tennullpunkte. Andererseits sind die Stellen der 
höchsten Biegebeanspruchung durch den größten 
Abstand des Kreises von der Wandungsmittellinie 
gekennzeichnet. 5 

. Die Längskraft, dividiert durch den Innendruck, 
erscheint als Länge des auf die Wandungsnormale 
projizierten, vom Biegezentrum zum Wandungs- 
punkt führenden Vektors. 

. Die Querkraft, dividiert durch den Innendruck, 
erscheint als Länge des gleichen, jedoch auf die 
Wandungstangente projizierten Vektors. 

Aus dieser Darstellung erklärt sich demnach auch 

in einfacher Weise die bekannte Tatsache, daß die 

Biegebeanspruchung um so größer wird, je mehr der 

Querschnitt: von der, Kreisform abweicht. 
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. Schließlich sei noch erwähnt, daß aus den Beziehun- 


en (8) und (9) auch für die Größe und Richtung der 
Fe kerondendchnittkereft einfacheGesetzmäßigkeiten 


- folgen. Wird der Abstand des Den 


unktes ' Biegezentrum mit r, und 
rischen r, und der Wandungsnormalen mitß bezeich- 
net, so folgt (Bild 3): 
2 = 14 008 (+ PB), y-yı = rrsin(a+P) (10), 
Demnach gehen Gleichung (8) und (9) über in 
Ip area al), 
g9g=—prrsnß.. vr. . (12). 
Für die Größe der resultierenden Schnittkraft folgt 
demnach 
ae DEE ee A LE): 


Ihre Richtung ergibt sich senkrecht zum Vektor r,. 
Unter Verwendung von R statt $ und Q wird zugleich 
_ das Kräftespiel anschaulicher. Wird wie in Bild 4 der 


Tl 


Homentennulur,f 


Bild 4. Kräftespiel am Winkelraum mit den Seiten o und 
7, bei Einführung der resultierenden Schnittkräfte R und 


rechte Schnitt in einen Momentenullpunkt gelegt, und 
die zugehörige resultierende Schnittkraft mit R, be- 
zeichnet, so folgt entsprechend Gleichung (13): 
R fi p 0 ar (14). 
Wird der von den Radien r, und o begrenzte Raum 
in die Gleichgewichtsbetrachtung einbezogen, so zeigt 
sich die translative Gleichgewichtsbedingung in an- 
schaulicher Weise von selbst erfüllt. Das Gleichgewicht 


' gegenüber einer Drehung entspricht der Bedingung 


M= 5 {n2 un 0?) AT NEE vl. (15), 


welche mit Gleichung (7) identisch ist. 


Die höchste Beanspruchung entspricht den extre- 
malen Werten des Radius r,, welche mit r% „., und 


r; min gekennzeichnet seien. Für die maximale Span- 
nung ergeben sich dann die Formeln: 

3 
Imax PD 5 (vr max 


bzw. 


3 : 1 } 
Omax =? E (e* Bin min) un Yen = (außen) (17). 


1 
e) + 7% "ymax| (innen) (16), 


Es bleibt noch anzugeben, wie die statisch Unbe- 
stimmten %,, 2, o am einfachsten ermittelt: werden. 
Hierzusei die Formänderungsenergie eines Quer- 
schnittstreifens (Breite in x-Richtung gleich 1) be- 
rechnet, Bei Vernachlässigung der Längskraftenergie 
folgt: 

U 
6 [M? 
na EZ) M% 
() 


ds 


+2 “ sr PerN: En “ | 
a en a 
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entsprechend der Integrationsvariablen | 
F2 r « N Be > 


DR Re 
er 


4 
8 4 
E K w' a  ı .-. [3 “ 
eingeführt. Die Anwendung des Satzes vom Extremum 
der Formänderungsenergie führt dann auf diedreiBe- 
ziehungen: s ER x 
04 04 9A 

” Er Fzciee Ra Fe We 7 (22), 

aus welchen 


4A 
oM oM °oM 
fr saan =, far anne [uno (23 
0 1) 


folgen. Mit Anwendung der Gleichung (7) ergibt sich 
hieraus: 2 


4- 4A 4A 
JM yar=0 [Mzdi=0, [Mdr=0 . (24). 
[) > 17 Ö 


Für die weitere Rechnung ist es praktisch, den Ur- 
sprung des Koordinatensystems y,z in den Schwer- 
punkt: des fiktiven Querschnittes mit der Wandstärke 
l/h? zu legen, d.h. derart, daß die Gleichungen 


erfüllt sind. Ferner mögen die Koordinatenachsen zu- 
Een a des Ersatzquerschnittes sein, so 
aß aue 


A 
[yzar=0 TR (26) 


gilt. Die Ausrechnung der Gleichungen (24) führt dann 
unmittelbar auf die gesuchten Beziehungen für die sta- 
tisch Unbestimmten. Werden zunächst Hauptträg- 
heitsmomente und Trägheitsradien des Ersatzquer- 
schnittes aus 


4 A 
3, = [ran 3,= [var y= = (m 


berechnet, so ergibt sich: 


4 
Da . 
= | Wer) ydR ... (28), 
(N) 
4A 


1? 
= 55, | WP+)zdi... .(), 

Ö h 
a (30). 
Der praktische Rechnungsgang besteht daher aus . 
folgenden Einzelschritten: 


1. Bestimmung des Schwerpunktes und der Haupt- 


achsen des Ersatzquerschnittes mit der Wand- 
stärke 1/h3; 
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Bild 5.. Beispiele. zur Biegebeanspruchung dünnwandiger 
Röhren unter Innendruck mit Angabe des Biegezentrums 
O,„ und des Radius. 


Gleichung (15), Ermittlung der größten und klein- 

‚sten Entfernung des Kreises von der Wandungs- 
mittellinie und Berechnung der größten Spannung 
aus Gleichung (16) bzw. (17); 

6. Zeichnen der normalen und tangentialen Projek- 
tionen des jeweiligen, vom Biegezentrum zur Wan- 
dung führenden Vektors, hieraus Bestimmung des 
Verlaufes der Längs- und Querkraft, entsprechend 
Gleichung (11) und (12). 

Einige Anwendungsbeispiele zeigt Bild 5. 
Dresden. - H. Neuber. 


Zur nichtlinearen Elastizitätstheorie.') 
Untersucht wurden die Beziehungen zwischen dem 
-Spannungstensor und Zerrungstensor (Deformations- 
tensor) in der nichtlinearen Elastizitätstheorie., Hierbei 
- wurden immer weitergehende einschränkende Annah- 
men gemacht. { 
AnnahmeTI!: Der Körper sei elastisch. Zu 
jedem Zerrungszustand gehört dann ein eindeutig be- 
stimmter Spannungszustand. Hieraus folgt,-daß jedem 
Zerrungszustand auch eine bestimmte eindeutige Zer- 
rungsenergie zugeordnet ist, die unabhängig vom Wege 
ı) Veranlassung zu dieser Mitteilung gab die kleine Mit- 
teilung von A. Philippidis, Z. angew. Math. Mech. 25—27 


(1947), 8. 31. BR 
2) Siehe z.B. A.E.H.Love, Lehrbuch der Elastizität 1907, 


8.69; die Komponente des. Tensors sind dort mit e,, und 
zy!2 bezeichnet. 


' gemeinen Zerrungstensor 


£ re. ne m‘ fs 

' Kreisprozeß Eı 
Kreisprozeß Energie gewinnen. PER 
Bei.großen partiellen Ableitungen der Komponenten 


u, v, w des Verschiebungsvektors erhält man den all- 
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Wir wollen am Anfange der Untersuchung diesen all- 
gemeinen Tensor benutzen. Für kleine partielle Ab- 


leistungen gehen die Komponenten des Tensors bei Ver- 


 nachlässigung höherer Glieder über in die Normalzer- 
rungen (Dehnungen) und Tangentialzerrungen (halbe 
Schiebungen). DE ET 


Aus dem Tensor (1) finden wir aus einer kubischen 
Gleichung die Hauptkomponente des Tensors &7, 8,» 
 &3g. Die Koeffizienten dieser Gleichung sind - . 


— (aıt &+ 3) =— (er +eyyt 
Ei1öga + Eggs t &stıı = Erz Tester 
ER | \ — (ey +... +...) 

——.&11 &93 CE ee 2 EL 


SET a8 Ei 

x Eyz Eyy®zx 77 ey)‘ 
Bei einer Drehung des Koordinatensystems bleiben 
die Koeffizienten invariant. Für unsere Untersuchung 


"nehmen wir diese Invarianten bzw. eine Kombination 


- derselben: 


Jıt&ıt &2t 83 =3emm Z 
Jg = (&11ı— Emm)? + (Eg2— Emm)’+ (Ess — Emm)? 
Se 
5 
J3 = E11 a2 Ey 


(3). 


YJi—3(&11&22 + &2&3 + &s &ı)] 


‚Sind die Hauptkomponenten 11, &33 und £,, des allge- 
meinen Zerrungstensors gefunden, so folgen daraus die 
Hauptzerrungen (Hauptdehnungen): 


4=V1+22,—1, 
oder 
221 =(l+.% —1, 28,=.. 


+ 
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DE 
‘Wir wollen nunmehr die mit den Zerrungen verbun- . 
dene Zerrungsenergie betrachten. Hierbei ist die Zer- 
rungsenergie A je Volumeneinheit des unverzerrten 
Körpers zu nehmen. Erteilt man dem.Körper für be- 
stimmte drei zueinander senkrechte Richtungen die 
Hauptzerrungen &,, & und &,, so ist bezogen auf die 
anfängliche Volumeneinheit eine bestimmte Energie 
vorhanden. Nimmt man dieselben Zerrungen, jedoch 
- für andere Richtungen, so wird allgemein auch die Zer- 
rungsenergie eine andere. Man denke z. B. an Holz, 
das einmal in Faserrichtung, ein anderes Mal senkrecht 
dazu gestreckt wird. Die Zerrungsenergie wird folglich 
eine Funktion sowohl der Hauptzerrungen &,, &, und g,, 
als auch der Richtungen dieser Zerrungen. 
Annahmell: Der Körper sei isotrop. Dann 
ist die Zerrungsenergie nur eine symmetrische Funktion 
der Hauptzerrungen e,, & und &,. Wir können hierfür 
aber auch die Hauptkomponenten des Zerrungstensors 


ee 


x &1, &y, und &,, und weiter die obigen drei Invarianten 


J1, J, und J, nach Gl. (3) nehmen: 
| A=A(J,J% J3) Cara Tr 7! = . (6). 


 Äus Ber Zerrungsenergie wollen wir die > ; 


nungen o,, o, und o, finden. Die Richtung derse 

fällt infolge der Isotropie mit der Richtung der Haupt- 
zerrungen &, & und &, zusammen. Wir betrachten ein 
Körperteilchen, das unverzerrt ein Würfel mit der 


. Kantenlänge 1 ist, wobei die Richtung der Kanten mit 


der Richtung der Hauptzerrungen zusammenfalle.. In- 
folge der Zerrungen werden die neuen Kantenlängen 
(1 + 8), (L+ &,) und (1 -+ &,). Wir beziehen jedoch 
die Spannungen auf die unverzerrten Anfangsflächen 
des Würfels. Dann wirkt in Richtung von e, die Kraft 


c, normal zur entsprechenden Fläche. Die Zerrungs- 


energie des Körperteilchens ist A. Vergrößern Y nun 
&, um de,, so vergrößert sich A um den Betrag 5, de,. 


1 
. Die Kräfte auf das Körperteilchen verrichten hierbei 


die Arbeit o,d e,. Hieraus folgt: 


9A 
Ran mV er ee (7) 


Hervorgehoben sei nochmals, daß die Spannungen und 
die Zerrungsenergie auf das unverzerrte Körperteilchen 
bezogen sind. Da A durch die Invarianten ausge- 
drückt ist, so folgt: 
942, DA Od, 2A 04, 
na Te ee 
Annahmelll: Die Normalzerrungen e,, &, und &, 


seien klein, so daß bei Vernachlässigung der höheren 
Glieder nach Gl. (4) gesetzt werden kann: 


Em am, mh. .n... (9). 
Dann ist 
J=4+8+&8=3en, 
I, = (&ı — Em)” + (&a — Em)” + (&s — Em)? 
2 (10). 


= (i HE Es — E28, — &g E53 — 8361) 


J=:8& 
Wir erhalten 


0A 0A 04 
ar an 2 (& — em) 57, + &2d7, (ıl). 


Hieraus 


1 
(m= 3 (1 +04 0,) 


ts atnata ar 
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AnnahmelV: Die Zerrungsenergie sei unabhän- 
gig von J,. In der linearen Rlastizitätstheorie 
wird die Zerrungsenergie A = (, J? + 0, J,. Für die 
nichtlineare Blastizitätstheorie wollen wir entsprechend 
annehmen, daß A unabhängig von J, ist, also 


Am AIET)  REN (13) 
Dann wird 
ß) A «Tr N Hr 
4 —-Om=2(e, — Ey)” Ä nn EN) 
9A (J,,J;) 
my, NEE (15) 


AnnahmeV: Der Ausdruck für A zerfalle in zwei 
Summanden, von denen der eine nur von /J,, der andere 
von J, abhängt: 


RUN EAN A EEE (16). 


Dann wird 


4 — Om —=2(E — Em) DIR Sehr) 
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Annahmen I bis V: 


Wir erhalten hiermit aus- 
Folgerung a: 


e v, 


om = Imlleı—em)t + (sem? + (sem (9) 


Folgerung b: 


£, 77 Em . 


Anstelle der Annahme Vkannmanauchfolgende 
Annahmen hen: SUR? 


mac 


Var Bee 


Aus jeder dieser Annahmen erhält man mit Hilfe der 


G.(15) bzw. (14) die Beziehung 
A=-Al)+4(d) 

und folglich dasselbe Ergebnis, wie aus der Annahme V. 
Es ist also gleichgültig, welche der drei Beziehungen als 
fünfte Annahme gewählt wird. 

Die Untersuchung zeigt, welche Annahmen in der 
nichtlinearen Elastizitätstheorie zu machen sind, um 
als Ergebnis Gl. (19) und GI. (20) zu erhalten. ; 


Schlewecke ü. Bockenem (20b). C. Weber. 


Zur Radizierung mit der Rechenmaschine. 

Hierüber ist ja schon öfter geschrieben worden, auch 
an dieser Stelle; meist handelte es sich um Quadrat- 
wurzeln. Die Frage, die heute hier beantwortet werden 
soll, ist vielmehr die, wie man mit einer Rechenma- 
schine am vorteilhaftesten Wurzeln mit beliebigen Ex- 
ponenten berechnen kann. Dabei ist zu unterscheiden 
zwischen der Verbesserung eines Näherungswertes und 
der direkten Berechnung einer Wurzel. 

Ist aus einer Tafel oder vom Rechenschieber ein 
Näherungswert a bekannt, so kann man mit der be- 
kannten Formel 


(r—1)0"+-.R 


na"! 


durch eine Division die Stellenzahl verdoppeln. Es soll 
nun gezeigt werden, daß es eine einfache und allgemeine 
Formel für quadratische Approximation gibt, die nicht 
mehr Rechenarbeit macht und dabei die dreifache 
Stellenzahl liefert. 

Dazu gehen wir aus von der binomischen Reihe 


N. 
YR= 


er As 1 n—l 2n—1)(n—1 
YyI+tz=1 + el Kb a; ee o—. 
2 |, .r—1 \, @r-1l)(n-1) 
ET a RE 
Nun ist bis auf Glieder dritter Ordnung 
ARSTER TE (> ) 
ES TE 
also bis auf Glieder 4. Ordnung 
n— RR pe de n—1[2n—1 n—1 
ae Nr: m 3 = ie 
2n 
und hieraus folgt für quadratische Annäherung mit 
R— a" 
= N 


die einfache und symmetrische Formel 


VR-a®R td R+m—ı1)a® 
Rl)R+M+Ha| : : 
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ai Der‘ Fehler ist kleiner as 5 
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Um die Stellenzahl der Maschine besser ausnutzen Ei x 


“_ können, wendet man für höhere. Genauigkeit die For- 


we ‚mel in folgender Form an 


nk 2 ee ee 
Ir VYR= 1 EEE u a EN . 
Se ern ni 
Ne Für m = 2 gehen a) und (2) über in R 
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Die Genauigkeit des neuen Näherungswert s a’ läßt. 


sich hier am einfachsten feststellen durch Bildung des 
BTRIE 1 R 

Quotienten R/a’; andererseits ergibt 2 (ar + = einen 
„verbesserten Näherungswert (Methode von Heron). 


Als Beispiel soll /2 berechnet werden. Mit dem Nä- 


herungswert 1,414 ergibt sich der auf 10 Stellen richtige 
Wert 1,41421 356236. E% : : CR RUE 


' Fürn=3 erhalten wir. 
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‘Als Beispiel soll v3 berechnet werden. Mit dem 2 
Näherungswert 1,442 ergibt sich der auf 10 Stellen 
‚ richtige Wert 1,44224 95703. 


Handelt es sich aber um die direkte Berechnung von 


Wurzeln, so wird man einfacher durch Anwendung der .. 


Formel 
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zum Ziel kommen; besonders für- große n wird sich 
dies immer empfehlen. Die hierbei notwendige Berech- 


nung des natürlichen Logarithmus durch inverse Inter- - 


polation der e-Funktion und ihre anschließende direkte 


Interpolation kann mit den vom Verfasser in seiner 
Arbeit „Praktische Interpolation höherer Ordnung‘ 
angegebenen. Formeln in der einfachsten Weise aus- 
geführt werden, so daß die Wurzelberechnung mit Ma- 
schine und Exponentialfunktion selbst für 10-stellige 
Genauigkeit als die rationellste bezeichnet werden darf. 


Göttingen. - . FR, Rubbert. 


_ BUCHBESPRECHUNGEN 


€. Weber, Festigkeitslehre. Wolfenbüttler 
Verlagsanstalt GmbH., Wolfenbüttel-Hannover 1947. 
Auf dem Gebiete. der Technischen Festigkeitslehre 
bringt das neue Buch insofern eine erfreuliche Bereiche- 
rung der Literatur, da es dem Verfasser gelungen ist, 
das an sich fachlich nüchterne Gebiet der Festigkeits- 


rechnung nicht nur in besonders anschaulicher Form . 


dem Leser näher zu bringen, sondern darüber hinaus 
auch in ein leicht verständliches, fast heiteres Gewand 
‘zu kleiden. Durch eine besonders klare und fast vor- 
"aussetzungslose Ausdruckstechnik, welche durch wür- 
zige und sprichwortartige Randbemerkungen noch 
unterstützt wird, führt der Verfasser auch den weniger 
vorgebildeten Leser unvefmerkt in jene Gefilde der 
technischen Festigkeitslehre hinein, welche sonst nur 
bei gewissen mathematischen Vorkenntnissen zugäng- 
lich sind. Auch hinsichtlich der begrifflichen Wort- 
prägung, der buchstabenmäßigen Bezeichnungsweise, 
sowie der Systematik in der Kräfteskizzierung be- 
schreitet der Verfasser in verschiedener Hinsicht neue 
Wege, die mancherlei Vorteile bieten und zu einer Re- 
form der seither in der technischen Festigkeitslehre 
üblichen Darstellungsart anregen. Neben der äußeren 
Form bringt das Büchlein auch inhaltlich manche neu- 
artigen Gedankengänge, so daß es nicht nur den Ler- 
nenden, sondern auch dem älteren Ingenieur zur Auf- 
frischung und Belebung seiner Kenntnisse Freude be- 
reiten wird. 


Dresden. H. Ne ob er. 


Dr. W. Lenz (o. ö. Prof. a. d.-Univ. Hamburg), Ein - 
führungsmathematik für Physiker. 
. (Büeher der Mathematik und Naturwissenschaften, 

herausgegeben von Dr. H. Poltz). 95 S. m. 45 Bildern, 
Wolfenbüttel-Hannover-1947, Wolfenbütteler Verlags- 
anstalt G.m.b. H. (Notdruck). Preis brosch. 8 RM. 


Die Notdrucke der Wolfenbütteler Verlagsanstalt 
wollen den heutigen großen Mangel an wissenschaft- 
licher Fachliteratur beheben. Sie sollen den Studieren- 
den die Möglichkeit zur Durcharbeitung desin den Vor- 
lesungen gebotenen Stoffes bieten und den bereits im 
Berufsleben Stehenden als Arbeitsunterlage dienen. So 
ist auch das vorliegende Heft, wie Verf. selbst sagt, ein 
Notbehelf. Es will den Studierenden auf möglichst 
kurzem Wege, wiederholend oder erstmalig erlernend, 
die mathematischen Kenntnisse zum Verständnis 
einer theoretisch-physikalischen Vorlesung vermitteln. 


"'Strengeoder allgemeine Beweise sind nichtgegeben. Das 


ist natürlich im allgemeinen nicht zu beanstanden, aber 
gerade bei den grundlegenden mathematischen Be- 


- griffen wünschte man doch eine eingehendere Behand- 


lung, die ohne wesentliche Vergrößerung des Umfanges“ 
möglich wäre. > = 

Manches, was sich sonst in einführenden Lehrbüchern 
findet, fehlt hier, so ist z. B. vom Differential, das auch 


, der Physiker gern gebraucht, nicht die Rede; dafür | 


werden gleich nach Einführung des Integrales Fourier- 


sche Reihen, Fouriersches Integral, danach erst Metho- 


den zur Berechnung von Integralen, unter ihnen die 
Sattelpunktsmethode, dann Integrationen im Kom- 
plexen und das Cauchysche Integral behandelt. Den 
Schluß bildet ein kurzer Überblick über Gestalt und 
Theorie der gewöhnlichen und partiellen Differential- 
gleichungen der Physik. 

Die Darstellung ist trotz der Knappheitleicht faßlich. 
Sie eignet sich besonders zur Wiederholung; dem An- 
fänger dürfte es an manchen Stellen kaum möglich sein, 
ihr mit Gewinn zu folgen. Auf jeden Fall scheint es mir 
aber nötig zu sein, daß nach dem Studium dieses Heftes 
die eingeführten Begriffe und Methoden in besonderen 
Vorlesungen sorgfältig besprochen werden, wobei dann 
auch auf die Methoden der angewandten Mathematik, 
die überhaupt nicht berührt werden, einzugehen wäre. 


Dresden. Fr. A. Willers. 


Raymond Claire Archibald, Mathematical 
table makers, portraits, paintings, busts, monu- 
ments, bio-bibliographical notes (New York. 1948, 
seripta mathematica), 82 Seiten, 17 cm breit, 24,5 c 
hoch, gebunden. 5 Se 


Der Herausgeber der Zeitschrift „Mathematical tab- 
les and other aids for computation“‘, Prof. Archi- 
bald von der Universität in Providence. (Rhode Is- 
land) hat eine Sammlung von Daten über 53 Berechner 
von Tafeln mit 20 schönen Bildnissen herausgegeben. 
Die 53 Mathematiker, Astronomen und mathemati- 
schen Physiker sind alphabetisch geordnet. Nach eini- 
gen biographischen Angaben folgen Mitteilungen über 
die vorhandenen Bildnisse und eine Aufzählung der 
wissenschaftlichen Arbeiten und veröffentlichten Zah- 
lentafeln. Die folgenden 12 deutschen Autoren werden 
behandelt: Anding, Bauschinger, E.E.H. 
Becker, Bessel, B.Cohn, Dase, Hoppe, 
Jacobi, Kepler, Lohse, Lommel, Pe- 
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ters. Das große Format des Buches ermöglichte es, 
N 


die Gesichter auf den ganzseitigen Bildnissen auf 
‘tem, glänzendem Papier deutlich erkennbar zu halten. 
- Die biographischen und sonstigen 


schichte der 


 EINGEGANGENE BÜCHER 


Bei der Schriftleitung sind folgende Bücher eingegangen (ausführliche Besprechung bleibt vorbehalten. 


D.. W. Weizel, (o. Prof. an der Universität Bonn), 
Einführungin die Physik. BandI: Me- 
chanik. (2. verb. Aufl.) Meyers kleine Handbücher 
Bd.42, 145 8. m. ır9 Abb., Leipzig 1947. Verlag: 
Bibliographisches Institut. Preis brosch. 2,90 M. 


Dipl.-Ing. K. Parnemann (Oberbaurat), DasRech- 
nen in der Elektrotechnik (Bücher der 
Technik. Herausgegeben von Dr.Ing.A. Kuhlenkamp) 
Teil I: Maßsysteme und Größenglei- 
chungen. 688. Preis brosch. €, —M. Teil H: 
Liniendiagramme — Zeigerdiagram- 
me. 67 S.m. 73 Abb. Preis brosch. 6,— M. Teil III: 

Die symbolische Rechnung. Ubungs- 
aufgaben. 60 8.m.41 Abb. Preis brosch. 4,80 M. 


NACHRICHTEN Er 


Max-Planck-Gesellschaft zur Förderung 
„. . der Wissenschaften. 

Am 26. Februar 1948 wurde in Göttingen unter der 
Billigung der britischen und der amerikanischen Militär- 
regierung die „Max-Planck-Gesellschaft zur Förderung 
der Wissenschaften‘‘ gegründet. Sie hat von den beiden 
Militärregierungen die Zusicherung erhalten, daß in ihr 
die wissenschaftliche Arbeit völlig frei und nur dem 

- Gesetz Nr. 25 unterworfen sein soll. Die erforderlichen 
Geldmittel sollen fürs erste von der Zwei-Zonenverwal- 
tung gewährt werden. 

Bis jetzt sind der Max-Planck-Gesellschaft ange- 
schlossen die Institute der früheren ‚‚Kaiser-Wilhelm- 
Gesellschaft zur Förderung der Wissenschaften‘ in der 
Britischen und Amerikanischen Zone sowie einige wei- 
tere Forschungsinstitute von wissenschaftlichem Rang. 
Der Beitritt von Instituten, auch solcher in anderen 
Zonen, ist gestattet, sofern die zuständigen Militär- 
regierungen ihre Erlaubnis hierzu erteilen und die Max- 
Planck-Gesellschaft zustimmt. Bis jetzt gehören der 
Max-Planck-Gesellschaft 34 Forschungsinstitute an. 

San% 


Nach dreijähriger Pause ist wieder im VDI-Verlag 
(jetzt Ratingen bei Düsseldorf) die Zeitschrift des Ver- 
eins Deutscher Ingenieure erschienen. Anfang März 
kam Band 90 Heft 1 heraus. 


Demnächst wird im Leibniz-Verlag (bisher R. Olden- 
bourg), München das ‚Allgemeine statistische Archiv“ 
wieder erscheinen, das unter anderem Arbeiten über 


Angaben beruhen 
offensichtlich auf einer jahrelangen gründlichen Quel- 
lenforschung. Dieses an ist nicht nur für die Ge- 
issenschaft von Bedeutung, sondern auch 
für die Menschenkunde. Die Bildnisse zeigen uns, daß 
Die besprochenen und angezeigten Bücher sind durch den 
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r 1947. Wolfenbütteler 
führungin die Ph sik. Band3 — 0, 
— Atomphysik — Wärme. 
Handbücher Bünd ER 155 8. mit 98 

Leipzig . Verlag Bibliographisches 
brosch. 2,90 RM. er 


Dr. Fr. Becker (o. Prof. an der Univ. 
führungin die Astronomie. 
Handbücher, Band4l. 148 8. mit 38 | 
2 Bibliographisches Institut. Preis 
2,90 RM. oz 
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Stuttgart. Am 13. Februar 1948 erlag der ord. 


Professor für Geodäsie und Direktor des Geodätischen Nu 
Institutes der TH. Stuttgart, Prof. Dr. phil. Edwin Ken 
Feyer im Alter von 60 Jahren während der Vor- Ag) 
lesung im Kreise seiner Studenten einem Herzschlag. e % 


Nach längerer Vorbereitung gründete das National Radar 
Bureau of Standards der USA. am 1. Juli 1947 die 
National Applied Mathematies Laboratories, 
unter der Leitung von J.H. Curtiss, die aus vier De‘ 
Instituten bestehen. irn | 

1. Das Institut für numerische an der 
Universität. des Staates Kalifornien in N k 
dieses Institut steht unter der Leitung von J. Todd, ne 
H.D. Huskey, O.Szasz und D.R. Hartree und "2 
wird im Laufe des Jahres 1948 seine Arbeiten auf- „A 
nehmen. - 

2. Das Rechenlaboratorium in New York oder” 
Washington unter A.N. Lowan, dem die Berechnung 
und Herausgabe von Funktionstafeln übertragen ist. 

3. Das Laboratorium für Statistik der Ingenieur- 
wissenschaften unter Leitung von C. Eisenhart, das 
die Anwendung moderner statistischer Methoden auf ( 
physikalische und technische Fragen untersuchen soll. 

4. Das Laboratorium für die Entwicklung maschi- 
neller Hilfsmittel unter Leitung von E.W. Cannon{ 
Es besteht aus einer mathematischen und einer elektro- 
technischen Abteilung und wird sich vor allem mit der 
Weiterentwicklung der automatischen Stellenrechen- 
maschinen beschäftigen. 

Der Etat der Laboratorien beträgt jährlich 532 000 


u 


mathematische und Versicherungs-Statistik bringen z' Dollar. Dazu kommt noch eine einmalige Summe von . 


wird. 


700 000 Dollar für die maschinelle Einrichtung. 
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